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MÉCANIQUE : CÉLESTE. 


NOTE 


SUR 


LES MÉTHODES DE WRONSKI", 


PAR M. YVON VILLARCEAU. 


Dans la Réforme des Mathématiques (t. 1, p. 163), Wronski distingue 
deux systèmes de Mécanique céleste, qu'il appelle science du désordre et 
science de l'ordre. Le premier est celui des analystes qui se sont occupés 
avant lui de la Mécanique céleste; l'autre est celui qu'il propose de substi- 
tuer à l'ancien, c'est-à-dire au système qui continue d'avoir cours dans la 
science. 

Ce géomètre a jugé utile de développer les deux systemes, qu'il caracté- 
rise par les équations différentielles auxquelles ils donnent lieu : le nouveau 
système est représenté par les équations (197), et l'ancien par les équa- 
tions (198). Il nous informe qu'il a fait l'application de l'un et de l'autre à la 
théorie de la Lune. Voici comment il s'exprime (p. 166, derniere ligne): 

« Nous avons ainsi employé les deux systémes pour offrir, d'une part, une 
immédiate vérification des résultats, par leur concordance dans ces deux 
systemes, et, de l'autre part, une preuve de l'extréme simplicité de la nou- 


(1) Lue au Bureau des Longitudes, dans sa séance du 23 mars 1881. Le Bureau a décidé, 
dans la méme séance, que cette Note et le Mémoire à l'appui seront insérés dans ses Annales; 
ladite Note a été publiée dans les Comptes rendus des séances de l'Académie des Sciences, 
t. XCII, séance du 4 avril 1881. 
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velle et véritable Mécanique céleste, par sa comparaison avec les inextri- 
cables complications de la fausse Mécanique céleste que l'on a poursuivie 
jusqu'à ce jour. » 

Un peu plus loin, il ajoute : 

« Tout ce que nous pouvons dire ici concernant la fausse théorie de la 
Science actuelle, c'est que l'orbite de la Lune, qui y résulte des prétendues 
perturbations causées par le Soleil, n'est nullement identique avec l'orbite 
variable que découvre la vraie théorie par l'influence téléologique du 
Soleil. » 

Deux faits ressortent de ces citations, l'un établissant l'accord des résul- 
tats définitifs des deux systèmes, l'autre faisant ressortir une grande diffé- 
rence dans les moyens d'obtenir ces résultats. Il est clair que l'expression 
de concordance, dont se sert Wronski, s'applique exclusivement aux expres- 
sions des coordonnées.de la Lune en fonction du temps. L'ancien et le nou- 
veau systeme seraient donc également exacts, et les résultats ne différe- 
raient que par la forme, circonstance qui néanmoins mérite toute l'attention 
des astronomes; puisque certaine méthode, se rapportant à l'ancien sys- 
teme, exigerait le calcul de deux mille termes, au moins, pour obtenir les 
trois coordonnées de la Lune à un instant donné. Nous devons exprimer ici 
le regret que le travail de Wronski, sur la théorie lunaire, soit resté à l'état 
de manuscrit, et de ne pouvoir opposer, à ce chiffre de deux mille termes, 
un autre chiffre de beaucoup inférieur et qui conduirait à des résultats 
bien plus pratiques. Toutefois, Wronski nous informe encore qu'il a fait, 
en cette circonstance, l'application de ses méthodes d'intégration : or les 
exemples qu'il a donnés d'autres applications des mémes méthodes peu- 
vent faire présumer des solutions avantageuses. 

Le second point que nous avons à examiner est celui qui concerne la 
grande différence des orbites dans les deux systèmes, différence qui sem- 
blerait tout d'abord inconciliable avec la concordance des résultats définitifs. 
Comme les géomètres de l'ancien système, Wronski fait usage de la mé- 
thode de variation des constantes arbitraires ou des éléments de l'orbite, 
mais il en fait un usage différent. Dans un systéme de trois équations diffé- 
rentielles du second ordre, on n'a, pour déterminer les variations des six 
constantes arbitraires, que trois conditions à remplir: celles de satisfaire 
aux trois équations différentielles du second ordre. Le probléme reste donc 
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indéterminé. Les géomètres ont depuis longtemps levé l'indétermination, 
en s'imposant la condition suivante : ils veulent que, par un choix conve- 
nable des constantes arbitraires, non seulement les expressions des coor- 
données aient la méme forme dans les deux cas de mouvement troublé et 
de mouvement non troublé, mais qu'il en soit de méme à l'égard des déri- 
vées premières des coordonnées ou des composantes des vitesses parallèles 
aux axes fixes. Ainsi se trouvent déterminées, pour les géomètres, les varia- 
tions différentielles des six constantes arbitraires. C'est là, assurément, 
une solution très correcte et très élégante; mais ne pourrait-on pas satis- 
faire aux mêmes conditions et d’une manière beaucoup plus simple, en 
choisissant convenablement les constantes? Tel est le but que s’est proposé 
Wronski, bien qu'il n'ait pas énoncé lesdites conditions. La préférence 
devra, ce nous semble, être accordée à celui des systèmes qui se prêtera le 
mieux aux intégrations à effectuer et qui exigera le moins de travail pour 


parvenir à un même degré de précision. 


Les méthodes de Wronski sont-elles exactes? On pourrait en douter, 
puisqu'elles n'ont donné lieu à aucune application de la part des astronomes 
ou des analystes. Il n'est pas nécessaire de rappeler les éloges donnés à ce 
géomètre par les commissaires de l'Académie des Sciences, lorsqu'ils ont eu 
à juger la haute valeur de la loi suprême des Mathématiques, pour justifier _ 
l'intérét que peuvent présenter les travaux ultérieurs du méme auteur. 

Nous ne pouvons croire à l'indifférence des astronomes à l'égard des pro- 
gres de la Mécanique céleste, et nous nous expliquons fort bien, les ayant 
éprouvées nous-mémes, les grandes difficultés qui s'offrent au lecteur le 
moins défavorablement prévenu contre les idées de Wronski. 

Après avoir, à plusieurs reprises, entrepris et abandonné l'étude des mé- 
thodes de Wronski, relatives à la Mécanique céleste, nous avons fini, cepen- 
dant, par nous rendre compte de ce qui peut avoir un véritable intérét, pour 
les astronomes-géomètres, dans les travaux publiés par Wronski sur cette 
matière. L'objet de la présente Note est d'exposer sommairement le résultat 
de nos recherches. 

Les deux méthodes, ancienne et nouvelle, ont cela de commun : l'auteur, 
au lieu de considérer les projections des forces perturbatrices sur des axes 
arbitraires et de directions fixes, s'attache à effectuer les projections sur 
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trois axes mobiles, dont l'un coincide avec le rayon vecteur, le deuxième 
situé dans le plan de l'orbite et perpendiculaire au précédent, le troisième 
perpendiculaire aux deux autres. 

L'objet principal du développement de l'ancienne méthode est unique- 
ment d'introduire, dans les équations différentielles du mouvement, les 
trois projections ou composantes des forces perturbatrices que nous venons 
d'indiquer. Nous n'essayerons pas de faire comprendre comment procède 
Wronski pour dégager les valeurs des différentielles des éléments de l'or- 
bite : nous avons dü renoncer nous-méme à le comprendre. Comme il 
importait surtout de vérifier l'exactitude des résultats, nous avons appliqué 
directement la méthode de la variation des constantes arbitraires. 

Voici ce que nous sommes parvenu à élucider à l'égard de l'ancienne 
méthode. Nous avons trouvé, en tenant un compte exact de la différence 
des notations, que les expressions différentielles du demi-grand axe, de 
l'excentricité et de la longitude du périhélie (') déduites par Wronski s'ac- 
cordent exactement avec les expressions en usage. Il n'est pas question de 
la longitude moyenne de l'époque, attendu que la longitude vraie, qu'elle 
sert à déterminer, est exprimée par une intégrale basée sur le principe des 
aires, que Wronski range avec les autres intégrales à traiter suivant ses 
méthodes. 

Quant aux éléments qui servent à déterminer la situation du plan va- 
riable de l'orbite, nous en trouvons, sous la marque (70) des Prolégoménes 
du Messianisme, des expressions différentielles où figure la moyenne arith- 
métique des vitesses de l'astre au périhélie et à l'aphélie. Au sujet de ces 
équations, Wronski fait cependant cette remarque |( p. 132): les expres- 
sions (70) « sont rigoureusement exactes, pourvu que, lorsqu'il s'agit de 
cette exactitude RAP on y substitue la vitesse vraie v à la place de la 
vitesse moyenne w... ». Or Je me suis assuré que les formules de Wronski, 
méme corrigées or ses indications, sont inexactes. Suivant nous, il fau- 
drait, pour en corriger l'erreur, multiplier les éléments différentiels de 


Wronski par le facteur 
Vi+ 2e COSY + e? 
1 + € COS 9 


(1) Les formules de Wronski contiennent la longitude de l'aphélie, au lieu de celle du 
périhélie. 
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où 9 désigne l'anomalie vraie et e l'excentricité. Ce facteur diffère très peu 
de l'unité, dans le cas des faibles excentricités, et il affecte la composante 
des forces perturbatrices suivant la normale au plan de l'orbite, en sorte que 
l'erreur commise, dans le cas de la théorie lunaire, peut se réduire à fort 
peu de chose; mais il y a là une erreur théorique dont l'origine nous parait 
étre dans la trop facile tendance de Wronski à s'appuyer sur des considéra- 
tions philosophi : nd cela n'est nullement nécessaire. Ces considéra- 
tions peuvent étre fort utiles pour asseoir solidement les bases d'une science 
qui n'est pas si absolument rationnelle qu'on aime à se le figurer; mais elles 
ne peuvent offrir ultérieurement d'autres avantages que de guider quel- 
quefois dans la recherche d'une solution difficile, et sous la condition de 
ne tolérer aucune dérogation aux principes admis. Comment donc Wronski 
a-t-il pu commettre l'erreur que nous signalons? Voici comment procède 
notre auteur. On lit (p. 291 des Prolégoménes du Messianisme) : 

« D'abord, pour établir le Canon astronomique, il suffit de considérer 
l'angle indéfiniment petit dp de la rotation qu'éprouve constamment, dans 
le temps indéfiniment petit dx (au lieu de dt), le plan de l'orbite sur le 
rayon vecteur r, par suite de l'influence principale de la susdite force per- 
turbatrice, dont l'action est perpendiculaire à ce plan... Or cet angle de 
rotation est... (suit la valeur de dg), et c'est là le principe ou la loi qui régit 
le Canon astronomique dont il s'agit. » 

On remarquera que Wronski ne démontre pas sa prétendue loi: il se 
borne à donner la formule qui la représente. C'est là le point de départ de 
l'erreur que nous signalons et que l'on corrigerait par l'application du fac- 
teur indiqué plus haut. Mais, nous le répétons, cette erreur ne doit avoir 
qu'une faible influence dans la théorie de la Lune. Ajoutons que les for- 
mules.qui règlent la Situation de l'orbite dans l'ancien et le nouveau sys- 
tème sont identiques. 


Passons actuellement à l'examen de la nouvelle méthode. Ici, comme dans 
l'autre, Wronski fait usage des composantes des forces perturbatrices sui- 
vant les trois axes mobiles. Pour établir les équations différentielles des 
trois éléments du mouvement, dans le plan de l'orbite, qu'il suffit de considé- 
rer, nous ne pouvons, par le motif déjà énoncé, exposer les démonstrations 
de Wronski : nous nous bornerons à indiquer les deux points qui nous sem- 
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blent caractériser la nouvelle méthode. Le premier est le résultat d'une 
simple remarque : l'action solaire, dirigée suivant le rayon vecteur, s'ajoute 
algébriquement avec la composante des forces perturbatrices suivant ce 
rayon; au lieu de rejeter cette composante, pour la ranger avec la seconde 
dans les termes de l'ordre des perturbations, Wronski la maintient jointe à 
l'action solaire. Le second point est relatif à la considération de la moyenne 
arithmétique w des vitesses aux deux extrémités du grand axe de l'orbite. 
Le produit de cette vitesse w et du demi-paramètre p est égal à la constante 
des aires. Telles sont p et w, les constantes variables que l'auteur de la 
nouvelle méthode substitue au demi-grand axe et à la constante des aires. 
Cela est évidemment permis. 

’artant donc de ces deux données, nous avons formé l'expression de la 
dérivée seconde du ravon vecteur par rapport au temps, et nous en avons 
déduit l'expression. rigoureusement exacte de la dérivée première, puis 
l'équation différentielle également exacte, de l'orbite, en coordonnées 
polaires. Admettant pour solution l'équation ordinaire de cette orbite, sous 
la condition d'y considérer les constantes comme variables, nous en avons 
déduit une équation de condition entre les variations du demi-paramètre p, 
de l'excentricité e et de la longitude w du périhélie. 

Cette condition étant supposée satisfaite, nous vérifions que le rayon 
vecteur et la longitude, ainsi que leurs dérivées premieres, conservent la 
méme forme, en fonction des éléments, que dans le mouvement elliptique. 

Une seconde équation entre les constantes est fournie par leur relation 
avec l'excentricité. Ces deux équations de condition permettent d'exprimer 
les différentielles de l'excentricité et de la longitude du périhélie, en fonction 
des différentielles de p et w; celles-ci s'expriment d'ailleurs au moyen des 
composantes des forces perturbatrices suivant les axes situés dans le plan 
de l'orbite, en sorte qu'il ne reste plus rien d'arbitraire : les différentielles 
des trois constantes p, e, 5 se trouvent ainsi entièrement déterminées. 

Nous avons comparé nos résultats avec ceux de Wronski, et constaté 
leur parfaite conformité. Comme dans l'autre méthode, notre géometre 
évite d'introduire la variation de la longitude moyenne de l'époque. Enfin, 
les éléments qui fixent la position du plan de l'orbite se déterminent par les 
mémes formules que dans l'ancienne méthode. 

Il suffit d'un simple coup d'oeil pour constater la simplicité remarquable 
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des expressions différentielles que fournit la nouvelle méthode comparée 


aux anciennes. 


Ici se bornent nos investigations. Elles établissent que, sauf une minime ; 
erreur facile à corriger, les formules de Wronski sont parfaitement exactes. | 


T —— n — 


Nous$erions heureux que cette constatation pùt déterminer quelque astro- 


nome-géomètre à poursuivre l'application des méthodes de Wronski à la 
théorie de la Lune. On y trouverait sans doute la solution des difficultés 
qui ont empêché les tentatives de Delaunay et de ses prédécesseurs d'abou- 
tir à un résultat complet, sous le double point de vue théorique et pratique. 

A ceux qui voudraient entreprendre un pareil travail, nous donnerions le 
conseil de consulter, s'il est possible, les manuscrits laissés par Wronski ('). 


A en juger par les travaux qu'il a publiés, on peut être certain de trouver | 
š rie £ g s | 
là une mine précieuse de développements analytiques, tres correctement | 


exécutés, On profiterait ainsi de l'expérience acquise par l'auteur de la nou- | 


velle méthode, et l'on se trouverait en présence de calculs faciles à vérifier | 
et à compléter, si l'auteur n'a pas entierement achevé son travail. 3 


(1) L'insertion de la présente Note aux Comptes rendus des séances de l'Académie des 
Sciences a été l'objet d'un article trés sympathique de la part du rédacteur du journal 
Wedrowieck (le Piéton), qui se publie à Varsovie. L'auteur y fait connaitre que la Société —— 
polonaise Atheneum se propose de provoquer la publication des manuscrits laissés par ` 
Wronski et qui sont la propriété de l'héritier du comte Dzialynski. 

Nous ne pouvons que faire des vœux pour qu'il soit procédé à cette publication dans le 
plus bref délai possible; nous estimons que les amis de la vérité et du progrés général des 
Sciences joindront leurs vœux aux nôtres. On saurait enfin, selon l'expression de l'auteur 
de l'article cité plus haut, « si Wronski n'a été qu'un grand charlatan, ou s'il était réelle- 
ment un homme de génie ». (Y; V.) 
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AVERTISSEMENT. 


La Note qui précède est la conclusion d'un travail de recherches concer- 
nant la substitution des composantes de la force perturbatrice suivant les 
axes mobiles avec le plan de l'orbite des corps célestes, aux composantes 
de la méme force suivant les axes fixes. 

Cette substitution a été mise en pratique par Wronski; mais l'obscurité 
dont ce géomètre a entouré l'exposé de ses méthodes n'ayant pas permis 
de le suivre, il fallait s'engager dans une voie conforme aux principes, uni- 
versellement admis, de la Mécanique céleste. 

Tel est l'objet des recherches que nous avons entreprises et que nous 
soumettons aujourd'hui à l'attention des astronomes. Ne pouvant prévoir, 
au début, quelles en seraient les conséquences, nous avons poursuivi ces 
recherches, sans nous astreindre à suivre un plan bien défini, mais avec la 
pensée d'utiliser les résultats successivement acquis, dans la recherche de 
résultats nouveaux. 

Ces quelques mots d'explication justifieront la forme sous laquelle nous 
présentons notre travail. 

Il est évident que nous aurions pu, après l'avoir terminé, en reprendre 
la rédaction, en vue des conclusions qui s'en dégagent; mais, outre l'incon- 
vénient de supprimer ainsi la véritable filiation des idées, il y aurait eu celui 
de retarder une publication dont l'urgence ne manquera pas de frapper les 
astronomes qui sont aujourd'hui arrétés par l'insuffisance des anciennes 
méthodes, lorsqu'ils ont à s'occuper du mouvement des nombreuses comètes 
et petites planètes dont les excentricités dépassent celle de la planète 
Mercure. 

On nous permettra donc de livrer notre travail à la publicité sans en 
modifier la forme primitive, et de nous attribuer la priorité d'une démons- 
tration intelligible des expressions différentielles sur lesquelles reposent 
les méthodes de Wronski concernant la Mécanique céleste. 


————— m GG c —— 
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NOUVELLE FORME 


DES 


ÉQUATIONS DIFFÉRENTIELLES 


pu 


MOUVEMENT DES PLANETES ET DES COMETES, 


PAR M. YVON VILLARCEAU. 


Le présent Mémoire a pour objet de mettre en évidence le mode d'in- 
fluence des forces perturbatrices, dans la variation des éléments des orbites : 
la forme que nous donnons aux équations différentielles est appropriée à 
l'emploi de méthodes d'intégration qu'il sera sans doute utile d'essayer, 
dans les cas où les méthodes en usage font défaut. Or ces cas sont très 
nombreux : la plupart des nouvelles petites planètes et toutes les cométes 
ont des orbites dont les excentricités excèdent les limites d'application des 
méthodes usuelles. Chacun sait d'ailleurs que la théorie de la Lune laisse 
beaucoup à désirer, et que Le Verrier a rencontré; dans celles de Jupiter 
et Saturne, des difficultés dont la solution analytique n'a pu acquérir le 
degré de précision habituel à ce grand astronome. 


Les géomètrés posent les équations du mouvement, en projetant, sur 
des axes fixes, les forces données et les accélérations. Par d'ingénieuses 
transformations analytiques, ils passent, des mouvements projetés sur les 
axes coordonnés, aux mouvements dans l'orbite elle-méme; et les résultats 


V. 2 
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qu'ils obtiennent sont trés simples et faciles à interpréter, lorsque l'astre 
considéré est soumis à la seule action solaire, les composantes des vitesses 
et des forces suivant les axes rectangulaires ayant disparu. Mais le probléme 
se complique extraordinairement, dés que l'on introduit les forces pertur- 
batrices : et il importe de remarquer que, si l'on s'affranchit encore de la 
considération des composantes de la vitesse suivant les axes fixes, en sub- 
stituant à la considération des composantes des forces perturbatrices sui- 
vant ces axes celle des dérivées de la fonction perturbatrice par rapport 
aux éléments, on ne forme pas les expressions de la variation des éléments 
en fonction des forces perturbatrices elles-mêmes. Il en est autrement, si 
l'on substitue aux composantes des forces perturbatrices suivant les axes 
fixes, les composantes des mémes forces suivant trois droites rectangu- 
laires, prises : deux dans le plan de l'aire élémentaire décrite dans le temps 
infiniment petit dt, parallèleme nt et perpendiculairement au rayon vecteur, 
et la troisième normalement à ce plan. 

Cette méthode a été signalée par Wronski, en 1842, dans les Prolégomènes 
du Messianisme et reproduite avec quelques développements, en 1847, par 
le méme géomètre, dans son Ouvrage intitulé Messianisme ou Réforme ab- 
solue du savoir humain. 

L'idée peut n'être pas nouvelle; mais la forme l'est assurément, en ce que 
celle adoptée par Wronski ne se rattache pas directement aux principes les 
mieux établis de la Mécanique, mais à des conceptions philosophiques, sur 
le róle des forces, qui sont inintelligibles pour les mécaniciens de notre 
époque. 

Si la forme offre de l'intérét, le fond doit être établi d'une manière inat- 
taquable : c'est ce que nous nous proposons de faire en ce moment. 

Parmi les divers procédés qui permettent d'introduire Îles composantes 
des forces suivant les axes mobiles que nous venons de désigner, nous n'en 
apercevons pas de plus simple, ni de plus direct, que l'emploi du théoréme 
de Coriolis, concernant les forces apparentes dans les mouvements relatifs. 
Nous regrettons seulement que cet auteur, préoccupé de la signification des 
forces apparentes, n'ait pas effectué le développement complet des termes 
auxquels elles donnent lieu. 

Cela nous met dans la nécessité de donner une démonstration directe des 
formules que nous aurons à emplover. 
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Soient æ, y, z les coordonnées d'un point matériel m, rapportées à un 
système d'axes rectangulaires fixes; æ,, y,, z, les coordonnées du méme 
point, rapportées à un système rectangulaire d'axes mobiles, ayant la méme 
origine fixe que l'autre système : nous donnons, dans le Tableau suivant, 
l'indication des lettres parlesquelles on représente ordinairement les cosinus 
des angles compris entre les directions des axes de l'un des systèmes et 
celles des axes de l'autre système. 


La fig. 1 établit la relation de position du plan des z,y, par rapport au 


Fig. 1. 


- 
Æ 


plan des æy, et celle des axes x et æ, par rapport à l'intersection N de ces 
deux plans. La première est définie par l'angle 6 des deux plans, les autres 
par les distances angulaires à et ọ de x et de x, au nœud N, le sens positif 
de ces trois angles étant celui qu'indique la figure. En admettant ces nota- 
tions, nous voulons utiliser les formules connues des mouvements de rota- 
tion : pour adapter finalement nos résultats aux usages des astronomes, dans 
la théorie du mouvement de translation des corps célestes, nous n'aurons 
besoin que d'y changer les signes des angles 0 et 4, en remplaçant en méme 
temps ces lettres par d'autres, ce qui évitera toute confusion. | 
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Le mouvement angulaire des axes mobiles étant entièrement arbitraire, 
nous assujettirons l'axe des x, à passer constamment par le point m. Nous 
voulons en outre que le mouvement effectif de ce point m ait lieu, à chaque 
instant, dans le plan des c, y,. Exprimons analytiquement cette condition : 
l'arc élémentaire ds étant situé dans le plan des z, y,, sa direction est per- 


pendiculaire à l'axe z, : on a donc la condition 
cos(ds, 3j) — 0; 
si l'on observe que les cosinus de l'angle de ds avec les trois axes fixes sont 


` dx fdy dz ; Å | 

respectivement —; =» =» et que (1) les cosinus des angles de z, avec les 
ds ÿds ds " 

mêmes axes sont c, c', c^, l'équation précédente donne 


c da + c'dy + c" da — o. 


Désignant le rayon vecteur par r, on aura suivant (1), attendu que la di- 


rection de x, est celle de r, 


(2) D'Or NET, FEU 
d’où l’on déduit 
(3) dx —adr+rda, dy—=a'dr+rda, da —a"dr + rda": 


mettant ces valeurs dans l'équation de condition, on aura 
(ac -- a! c' -- a" c") dr + r(cda + c' da! + c" da") — 0; 


or, le coefficient de dr est nul en vertu des relations qui lient les neuf cosi- 


nus : nous avons donc la condition 
(4) c da + c! da! + c" da" — o. 
On trouvera ( Mécanique de Poisson, 2° édition, t. II, p. 123) que la va- 


leur du premier membre de cette équation est — gdt, q désignant la com- 
posante de la vitesse de rotation du système d'axes mobiles autour de l'axe 


des y,. Notre équation de condition se réduit ainsi à 
(5) qq = 0. 


Nous allons actuellement écrire les équations du mouvement d'un corps 
céleste en y introduisant directement, à la place des.composantes des forces 
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suivant les axes fixes, les sommes X,, Y,, Z, des projections des composantes 
suivant les axes mobiles. En convenant que ces composantes se rapportent 
à l'unité de masse, nous aurons les trois équations 


d? x VV 

—— — 0 Xx.2-5 Y HP 
dg gd it € £y 
d? y V r r 
dB =“ A+ bV Y+ cL, 
d z E : 

— =X b Y,-2- c" Zi 
dt $t ` i 


Multiplions ces équations par a, a’, a” et ajoutons, nous obtiendrons, en 
vertu des relations entre les cosinus, la premiere des équations suivantes; 
opérant de méme au moyen des autres cosinus, b, b', b^, c, c', c", il viendra 


dx 3H dy E dz 
; ü 3 


OA —— r — ÓÀÀ e$ X 
dt? de de 2 
d*a d? y dæ sz s 
v xen c r INS 
dt? dte? dt? 
s diaz F d? y inier [n AE 
de de dit RN 
Or nous déduisons des équations (3) 
gue. Á [æi i da dr > dia 
di 7° ga dtu. ge 
æ y ái dir Pac da dr T d? q' 
PES k——- = D res, es Ji aenn a 
at? dte "t dt? 
U TA 4" dr AH da" dr De d? a" 
d? ^ dg dt dt dt? 


Substituant ces valeurs dans les équations précédentes, et observant que 


(6) a da + a! da! + a" da" — o, 
on aura 
da d? d? d? a v 
dt? Gi (a de dt? de Kas 
FA pda , , da , da” 1, da , Pa ,u * 
De hak c cris ei reb sp Y, 
al NS m NIA „ da" („Pa „ Pa , c s )- 7 
\ ail "dt aw T CEN CU nes Wu 
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Il nous reste à exprimer les (  ) qui figurent ici, en fonction des compo- 


santes p, q, r dela vitesse derotation (nous surmontons les lettres ret p d'un 
trait, pour éviter toute confusion avec le rayon vecteur r et avec le demi- 
paramètre p). En premier lieu, l’on déduit de (6) 


ada + a' d'a! + a" d* a! — — (da? + da? + da"); 


or les formules du mouvement de rotation (p. 135 du méme Volume de 
Poisson) donnent aisément 


da? + da"? + da'? — (r*+ q?) dt?, 


résultat qui se réduit ici à ? df, à cause de (5) : la première de nos équa- 
tions (7) devient en conséquence 


2 2 UE A 
(8) Br Xe 


D'un autre côté, les mêmes formules de Poisson donnent, à cause de g = o, 


b da + V da! + b da" =r dt, 
cda- c! da + c" da" — o; 


on en déduit, par la différentiation, 


b d*a + V/ da! + b' d? a" — — (dadb + da! dU' + da" db") + dr dt, 
c d* a 4- c! da + c" d* a" — — (da dc + da! dc' + da" dc"); 


or les équations (8) de Poisson (méme page) donnent 


da db + da! db! + da" db" = — pq dt, 
da de + da! dc! + da" dc" = — rpdt*. 


Au moyen de ces relations, et à cause de q — 0, nos deuxième et troisième 
équations (7) deviennent 


| Ao dry 
rrp = Li. T 


Rappelons actuellement les relations qui lient les composantes de la vi- 
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tesse de rotation avec les variations des angles 0, o, Ý : ces relations [équa- 
tion (7) de Poisson, p. 134] deviennent ici, à cause de g = o, 


D= sin un St S Lu 
IURE OS Pg? 


! 


o = coso sinô ds H sii ay 
= 2 dt In © dt? 


La combinaison des deux premieres donne 


d So ^ 
di^ P959 
(10) 


| sin d p sino; 
\ 


quant à la troisième, observant que r est la vitesse de rotation autour de 
l'axe des z, perpendiculaire au plan de laire élémentaire décrite par le 
rayon vecteur dans l'instant dt, cette vitesse r est la vitesse angulaire du 
rayon vecteur lui-même. Donc, si nous désignons par d6 l'angle de deux 
rayons vecteurs consécutifs, nous pourrons écrire 


= dde 
dt? 
et nous aurons 
dé de dy 
= -——cos0—. 
(11) dt dt dt 


Nous allons remplacer r par 7- = ; dans nos équations (8) et (9) et distinguer, 


dins la composante X,, les cos provenant de l'action du Soleil et les 
forces perturbatrices qui proviennent de l'action des planètes : la première 
de ces forces, agissant précisément dans la direction opposée au rayon vec- 


teur, sera — f E, en désignant par f le coefficient de l'attraction universelle 


et faisant 
(12) p=M+m, 


somme des masses du Soleil et de la planète troublée; quant aux autres 
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forces, nous en désignerons par X, la somme des composantes suivant r 
ou æ,. L'action du Soleil étant perpendiculaire aux axes y, et z, ne donne 
pas de composantes suivant ces axes : les trois lettres X,, Y,, Z, désigneront 
dès lors uniquement les composantes des forces perturbatrices. En consé- 


quence, nous aurons 


dir dp? fv. x 
aulas SC M 
jS Á dr dd do =Y : 
(13) 2] di 5 be. 1* 
TI Doe 
pr^ ~ Z 


Si nous multiplions la deuxième de ces équations par r, son premier 


membre deviendra 


, d^ 

a 

d.r? d ,d9 lt 

R f 
di di dË dii. 
et l'on aura 
, dẹ : 
l. »; ET paa Yit: 


Soit 


(14) | k= f Yiedr, 
nous aurons, en intégrant l'équation précédente, 
(15) r? — =k; 


c'est l'équation des aires. 
$ Pt S > k se cdd 
Portant dans la première et la troisième des équations (13) la valeur de se 
] € 


qu'on déduit de la précédente, et transposant, il viendra 


bos À Æ 1 
( 16) Øk æd e AE EE 
s. dt? p? D 
(17) ps i 


Au moyen de cette derniere, les équations (10), qui reglent le mouvement 
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mouvement du plan de l'orbite, seront 


1 SR Wm r 
| D ap 77 COS ©, 
RS / Z 

v AY T AR 
sung = - Sin o. 
sind 77 =+ -y sing 


Les équations (14), (15), (16) et (18), jointes à (11), sont celles qu'il reste 
à résoudre. 

L'équation (14) nous montre que les variations de k, ou du double de l'aire 
décrite par le rayon vecteur par unité de temps, sont entierement dues à la 
composante Y, des forces perturbatrices. Si l'on observe que la variation 
de £ ne produit, dans les expressions (18), que des variations du second 
ordre, on reconnaitra que les changements de l'inclinaison et de la longi- 
tude du nœud sont principalement dus à l'action de la composante Z, des 
forces perturbatrices suivant la normale au plan de l'orbite. 

L'emploi ultérieur de cette composante nous conduira à faire usage des 
valeurs des deux cosinus b et b en fonction des angles w, Ø et 4. Pour ne 
point avoir à y revenir, nous reproduirons ici ces valeurs, que l'on trouve j| 
dans la Mécanique de Poisson (t. II, p. 64) : 


(19) b' = cos0 cosy cosy + sin y sinc, 
(20) b = cos0 sin y cose — cos sin y. 


Intégration de l'équation (16), en négligeant les forces perturbatrices, 
c'est-à-dire en y supposant k constant et négligeant X,. 


Bien que cette intégrale soit trés connue, il ne sera pas sans intérét de 
l'effectuer; nous reconnaitrons plus clairement l'influence des composantes 
des forces perturbatrices sur les éléments que nous avons à déterminer. 

Multiplions par dr l'équation (16), réduite comme il vient d'étre dit, et 
intégrons ; nous aurons, en désignant la constante par — i E 


ou 


(21) 
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Avant d'aller plus loin, nous ferons remarquer que l'équation (15) donne. 


r? dẹ? DE KT 
nn d?  r?? 


si donc nous désignons par ds l'élément de l'orbite décrit dans le temps dt, 
nous aurons, en ajoutant cette équation et la précédente, | 
(2 

(28) F 
relation trés simple entre le carré de la vitesse ^ le rayon vecteur r et la 
constante U, relation sur laquelle nous reviendrons. 

Au lieu d'intégrer immédiatement l'équation (21), il sera plus simple 
d'obtenir d'abord l'équation polaire de l'orbite; ce qui se fera en éliminant 
dt? entre (21) et (22). Divisant donc ces deux équations membre à membre, 


il viendra 
2 
1 dr afe. — Per U r? 
r? db? k? Kt 
or, si l'on fait 
k 
24 FI 
(24) 2 
d'oü 
dr ita do 
Pr E 


notre équation deviendra 


de?  2fy1 U 
e? dé* - KK on SEE p? 
ou 
dè fiy foy 
m de Me (e- 


On intégrera cette équation en faisant 


ed 
k t 
(26) cost m AERE oU d'où e md sin £ = A P EE ie 
2132 addo dp 2 42 
fv Pr [^y T 
dod hb s 


La substitution de ces valeurs dans l'équation précédente donne, sup- 
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pression faite des facteurs communs, 


d’où 
C = + (P ksa v), 


en désignant par w la constante d'intégration. Mettant cette valeur dans 
l'équation (26), on en déduit 


dins: 
P Zb e vbt — U cos ($ — x) 


puis, en vertu de (24), 
k? 


RES, 


mum EP 
k 
+y gea Ueos(o — m 


équation qui est celle d'une section conique. 
Si l'on désigne par a le demi-grand axe et par e l'excentricité, on aura 


k? k? 


on) M =y REY AN ES . 
(27) a (1 UM E e!—1 Fu Vi 
d'où 

(28) aA 


et l'équation polaire de l'orbite prendra la forme ordinaire 


~ a(1— e?) 


(29) ^ r--ecos(9—w) 


La valeur de e? montre que U sera positif dans le cas de l'ellipse, nul 
dans le cas de la parabole, et négatif dans celui de l'hyperbole. 
Si l'on met dans (23) la valeur de U que fournit l'équation (28), on aura 


l'équation des forces vives 
À ds? EE. r^ 
(30) | =se- a) 


Revenons à notre objet, l'intégration de l'équation (21); multiplions 
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y ee R ; Ur? 

immédiatement l'équation (21) par fn nous aurons 

(31) Urt Uk pra 
| pe de pq) 


ou, en avant égard aux relations (27) et (28), 


r? dr? : rS? 
= € — | 1 — 3 . 


afy de? a 
Posons 
< r 
(32) I — — = e COS, 
a 
d'oü 
(39) "—1—ecosu, ét 4 e esinu su 
5r ib M ` RSR ne Dee 
a ú a dt dt? 
dr 


nous aurons, en substituant les valeurs de r et z; que fournissent ces rela- 


tions dans l'équation précédente, suppression faite des facteurs communs, 


æ du? | 
— (1 — e cosu) R=! 
fv dt 
Faisons, suivant l'usage, 
(34) ant zm yu 


et portons cette valeur dans l'équation précédente; on en tirera, après avoir 
pris la racine carrée des deux membres, 


(35) ndt = (1— ecosu) du: 


si l'on n'a pas mis de double signe, c'est afin que l'auxiliaire u et le temps t£ 
croissent simultanément. Intégrant et désignant par 7 la constante, il vient 
(36) n(t—*)-u--esinu. 

Pour éviter de certaines difficultés dans le mouvement troublé, on sub- 
stitue à cette formule la suivante 
(37) fndt 4-  — 5 — u — esinu, 


oü la somme des deux premiers termes désigne la longitude moyenne à 
l'époque 7. Par la premiere formule (33), on voit que le minimum de r a 
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lieu quand-cosu est égal à l'unité ou sin = o: il s'ensuit que la constante + 
désigne l'époque du passage au périhélie, ce qui justifie l'introduction de la 
longitude z du périhélie dans le premier membre de (37). 

De la combinaison des relations (29) et (33) on déduit un grand nombre 
de formules, dont nous rapporterons ici celles que nous aurons à utiliser, 
savoir : | 
COS U — e e + cos( — v) 


38 COS(P — 5) = cosu = 
i ( dc par 1+ecos(® — v) 


| Vr sinz (9 — v) — Va(14- e) sin > u, 
(39) 


| Vr cos = (4 — v) = ya(1 — e) cos= u, 


(ho) (n sin ($ — 5) — avy1— e sinu, 
| rcos(ó — w) —a(cosu — e), 
V1— e sin (£ — s) cosu — cos(® — m) sinu — e sinu, 


Vi 


V 1— e cos (€ — w) sinu — sin (# — v) (cosu — e) = o, 


e? sin (9 — v) sin u + COS(P — w) (cosu — e) = 1 — e cosu, 


1— e? cos (P — m) cosu + sin (P — m) sinu Vi. e, 


e : sÍ 
Fe (P—w)sinu + COS(P — m) = COS, 
1 — e? 
sinu 


e , ; 
—=— c0s(9 — w) sinu — sin(* — wv) — — — 


1— e? Vie 


‘ 


Nous ne poursuivrons pas davantage le développement de formules qui 
sont suffisamment connues; nous allons aborder la formation des expres- 
sions différentielles des éléments de l'orbite. 


VARIATIONS DES ÉLÉMENTS DE L'ORBITE (!). 


Nous ferons usage de la théorie de la variation des constantes arbitraires. 
Rappelons, en quelques mots, en quoi elle consiste. On considére les expres- 


(1) Les développements suivants se rapportent à ce que Wronski appelle la science du 
désordre. 
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sions des coordonnées, obtenues dans le cas du mouvement non troublé, 
comme applicables au cas des perturbations, moyennant l'emploi de valeurs 
convenables des constantes, qui deviennent dés lors de véritables variables. 
On assujettit ces nouvelles variables à satisfaire aux équations différentielles 
du mouvement troublé, et l'on obtient des équations en nombre insuffisant 
pour déterminer les valeurs de tous les éléments. Cette circonstance 
permet de poser de nouvelles conditions, qui consistent en ce que les expres- 
sions des dérivées des coordonnées, par rapport au temps, restent les 
mémes que dans le mouvement non troublé. Ces conditions sont faciles à 
remplir et suffisent à la détermination des valeurs variables des éléments. 
Nous n'avons à traiter ainsi que les éléments entrant dans l'expression 
des coordonnées ® et r, et qui sont s, a, e, w, attendu que les équations (11) 
et (18), en supposant qu'on en élimine la variable ©, permettront d'obtenir 
directement les valeurs de 0 et y. 


dr dde 
Ceci posé, nous allons former les valeurs des dérivées zi et > en faisant 


à la fois varier le temps et les constantes. 
Différentions les équations (40) : nous aurons 


sin(9 — w) dr + r cos ($ — w) d ($ — v) 


sinude, 


— avy1— e*cosu du + V1 — sin u da — 
` 2 
yı—e 


cos( — w) dr — rsin ($ — w) d ($ — v) 


= — a sinu du + (cosu — e) da — ade; 


on en déduit d'abord 


dr=  a|y1— e sin(& — w) cosu — cos(® — w) sin u | du 


+ |V1— esin(t — v) sinu + cos (P — v) (cosik- e)] da 


[ 
- y£ =. sin (P — w) sin uw + cos(* — o) |a, 


rd(9 —9)—  a|yV1— è cos(? — w) cosu + sin (£ — w) sinu] du 


+ [V1 — ei cos (P — w) sinu — sin (P — w) (cosu — e)] da 


cos (P — w) sinu — sin ($ —=»)| de. 


caes 
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En ayant égard aux relations (41), (42) et (43), ces expressions de- 
viennent 


diss LA C EH 2 de 
dr — “e sinu 77 + (1— ecosu) 75 fear d 


L^ dece dt 
rd(6—w) _- ac asinu de 
dt cae i a dt 


or, la différentiation complète de l'équation (37) donne 


du : 
—— = (1—eCoSu) — —sinu 
dt 


V/v 


d'où, en remplaçant 1 — ecosu par sa valeur = Letn par sa valeur =- Sarr tirée 


e (34), 


a 


du. y V fv. 1 de — dew dde 
= m a) qo simu o. 


Mettant cette valeur dans les équations précédentes, on aura 


5 A Fa tSimu æesinu eser) ae i SS: db ee 
Va dotem dt HV amy 0su) 7; e 
do  .-Vya(i—e) a de dọ ay T:— M ` asinu \ de 
MU se WAT 72 =a r? yı € (z— di e HET 7 Sin u + Mise — 
Si Pon élimine de ces formules l'anomalie excentrique u, elles deviennent 
dr /— wvfe ae v (á T) 
— = — D) + —— SIn ($ — o) | —— 
dt Va(1— e) ) ) Vi— e ( dt^ dt 
r da 
— a COS ( P — =) T+ 2t 
(44) 


EU jð a. — — de a? ; V do 
=V fe -Si-e$-(iI-Aac9)m 


sers. 1 de 
r — 


+ sin(9 — z) 


Ces équations fourniraient les expressions des dérivées de r et ® par 
rapport aux constantes e, æ, id a; mais tel n'est pas notre objet. 


db 
Les parties des valeurs de 47 ; etm j? qui répondent au mouvement non 
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troublé, sont les premiers termes des deux expressions précédentes; 
comme ces expressions doivent rester les mémes dans le mouvement trou- 
blé, nous devons égaler à zéro l'ensemble des termes qui viennent à la 
suite. On a ainsi 


= = dece — ©), 
A^ l Va(1— e 
(45) 
VR 
T= ng 


puis, en multipliant les termes restants de la premiere équation (44) par 


/1 F t 
A » et ceux de la seconde par ——— + 
a? V1 — e 
de la —— y. yp d 
esin(d— w) ^ — esin( —w) 77 — =æ cos(b — m) À + 2 1— e j 
Ee dt 
Y de r? dw r? sin(ó —w)[a(t—^e? L6 
(46) ata -Um "AN i 
at a wie at a V1— e! b. 7 e dt 


DN de £ : 
Eliminons d: entre ces deux equations, nous aurons 


ca n k e sin? (p — m) [E72 


T2 (1—e)yi—e r 


, 
d Pole dde TUE 


ou 
r ?sin(P — w) edw r? de r ,———1da 
-+ r [2e + (1 + e?) cos(b — > e - — 
a? V1 —e di ayi— e al C ) Nes T. ^t à V "ad 
d'oü 
1 da n I' ue 
_ EE sin (D — Et. 2 ex 
a dt  a(1— 5 | ( UE p +[2e+(1+ e) cos( s 1— e? dt | 
équation qui peut s'écrire 
(47) 1 da — sin(P—75) edw i 2e+{(1+e )cos(d —m) 1 de 
^ a dt `` 1--ecos(? —w) dt I + ecos(^ — s) r— ei dt 
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de 
De l'équation (46) on tire cette valeur de 7; 


; E 

di (1—e)? dw EAE =). de 

add" ' HR [.— € - a (V7 Ti recos(®— ©) mur 

Nous transformerons cette expression, quand nous aurons formé les valeurs 
dw de 


de e 7; et 7, 


Nous allons maintenant effectuer le calcul de la deuxième dérivée de r. 
Par la première équation (45), nous aurons 


dr Vfu d^ dw 
di = kal C a ar) 


V fv. I UE a 
—————— n A æg] 


sin(4 —«) Z + 3i 


va(r— e?) a? (1— e)? 


+ 


Le premier terme de cette expression peut s'écrire, en vertu de (29), 


| Vfeaü—e) — VÆ lE Y: A 


r Va(1— e!) di} 


de r Q , : 4 
Les termes en di reunis Se réduisent a 


mec a c de. 
Fee à 


va(1— e?) 
da 
le terme en T est 
I Vv. esin pesin( — w) 


) 3 de 
2 icu a dt’ 


ou, en y mettant la valeur (47), 


i NA EL esin?(?—w) edw 
Pa MER 2e + (1 + e?) cos( — m) » ı de 
2 Va(1— e) 1 + ecos (9 — x) Go A pd Dee —e dt 
4 
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Au moyen de ces valeurs, la dérivée dia devient d'abord 


= [VIR ai øði Vie. ap: 
E T 


r Va(1 — e) | dt 
4 Vf». cos( — s) "m kc AM. _ esin'($—w). T 
Va(1— e) g Va(1— e) I + ecos( P — w) 
V/v 1 2e-r- (1+ e?) cos(&. — v) Pv de. 
pe 2 I + ecos ( — m) e | since ra = dt 


de 
Or, si nous introduisons, dans le premier terme, la valeur de 7; k de T° et 


si nous remplacons en E V fp.a(1— — e*) par sa valeur £, suivant (27), 


le premier terme ded ds " deviendra 
KL: 
nB. pM 


edw de 
les termes en prs et + sont réductibles à 


1! AGE e + cos(P— m) 
1 VIS 2 + ecos( — 


2 DS v ee b — 
2 Va(1— e) ! + ecos(? — sj sin 25 di’ 


Sks cos (a — w) | us 


on a donc finalement 


Pr fW Wes V/v “e + cos($ — v) BEEN 
| # mp 2Va(1— e) 1 e C08 (P —T) 
(90) 
pr. Wi" 2 + e cos(P — v) 


ba S 
2 Va(1 — e) I + e cos ($ — v) dt 


En comparant cette équation à l'équation différentielle proposée (16), on 
voit que la condition de leur concordance est 


| 1 R fetes o que 


i 2 va(1— e?) I + ecos( P — w) 
ss | Va (9 — d 
I fv 2 + e cos r^ oes Ge REY 
| A A, dan Euro E (sr) sin (9 —w) 4 Se 
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Il reste à satisfaire à là seconde équation (14): or cette équation revient à 


1 da 
ou, en mettant pour = 7; sa valeur (47), 


2dk 3 smiS uw) eds |  cos(b—w) de 
k dt ^ 1--ecos(P —w) dt ' 1-- ecos(9 w) dt 


Au moyen de la valeur de ia on déduit de cette relation 


(5a) __ Sin(9 — v) e cos (P — x) de 2v, 
1--ecos(P —w) dt 1-r-ecos(P—cw)dt ` k ` 


d 
Les équations (51) et (52) vont nous donner les valeurs de æ et D 


i à " 1 da ; 

qui, étant mises dans (47) et (48), feront connaitre - apt ur 
Pour effectuer ces diverses déterminations, nous remplacerons dans (51) 

la quantité mr ar — k, suivant la relation (27), et ferons disparaltre 
q RO PA S dr P 


le dénominateur 2 [1 + ecos(® — &)]; nous aurons ainsi 


— [e + 2008 (P — v) + ecos? ($ — w)] - 
T [2 + e cos ($ — v)] sin (P — w) z = 2[(1 + e cos ($ — w)] kA 


En faisant disparaitre le dénominateur commun du premier membre 
de (52), on aura de méme 


sin(? — w) 42 + cos(& — m) Z Uu” Fl! + ecos(6 — «)] Yr a Zt» 


équations où il ne faut pas perdre de vue que £ représente l'intégrale (14). 


Nous allons tirer de ces équations les valeurs de wis et $ 


dt dé 
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Soient, pour abréger, 


(53) M-—[2-r-ecos($ — v)]sin (9 — x), N = e-- 2cos (P — Tv) + ecos? (t — w), 


(54) ALT, B—1Yr; 


ces équations deviendront 


de e dw 
Ma ua 


I] 


2 A[1+ ecos(P — w)], 
cos (P — v) = + sin (P — w) r = 2 B[1 + ecos(? — w)]: 


on tire de là 


—" 


[Msin(& — wv) + Ncos(b — v) = &[MB — A cos(% — w)][1 + ecós (^ — v)], 


[Msin(9 — s) + N cos(4 — «)] 5 = 2 [NB + A sin ($ —5)][14- ecos( — »)]; 


— 


or, le facteur commun de nos deux inconnues est égal à 2 [ 1 + ecos( — cJ]; 


il vient donc simplement 


aa = MB — A cos (P — w), 
dt 
(55) 
zd = NB + Asin ($ — x) 
dl T j 


Substituons ces valeurs dans expression (47) multipliée par [ t +ecos(d — a)]; 


nous aurons d'abord 


2e + (1+ e?) cos (P — =] 


ı da ù 
[1+ ecos($ — 9)] 2 77 = B [ Msin(e — v) + N Eme 


+ À [nne eom — cos (b — d sin ($ — w). 


Si l'on réduit au méme dénominateur les deux termes du facteur de B, 
on trouve, pour le coefficient de B, 


2 + 6ecos ($ — w) + 6e? cos? ($ — v) + 26? cos? (P — v) 2 
aisi duka A eese fra nes [1 + ecos(b — w)J*. 
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Quant au coefficient de A, il se réduit immédiatement à 


2e : 
yr [1 + ecos(9 — w)] sin (# — s); 


on a donc 


ae a. 
: sin (P — w) A. 


da 2 


I 
a 
Si l'on remplace l'un des facteurs [1 + ecos(® — w)] de B par sa valeur 
2 
star que l'on mette les valeurs (54) de B et de A, on aura 


fer 


(56) 


MÀ 2k 1 
Accra AU ee CAP cvs 


Transportons encore les valeurs (55) dans l'expression (48), multipliée 
par [1 + ecos(® — z)]?, nous aurons 


[1 + e cos (P — s)P? - 


e 


P ne el e Liens soni] u) 


3 
2 


ET [1 n (1 — ety: costa à) 


+Vi—e [2 + e008(8 — w)]sini (9 -t |. 


Il se présente ici des réductions. En vertu de (53), le premier terme est 
égal à 
3 
2 


BM Mix aea cui ewm) 


= Ex | [1 + ecos($ — s)? — V1— e (Ne + 1— et)j. 
Or, Ne + 1 — e? est égal à [1 + ecos(® — ø )]?; le premier terme est donc 


BM BMe 1, 
TT [1+ ecos(® — s)? (1— /1— e) = FY nme [1 + ecos(9* — »)]'; 
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le second terme se transforme en 
À 
— zit + e cos ( — æ )]?cos(P — v) 
+V1— e | [2e-- e*cos($ — s)][1— cos? ( — s)]— (1— e?) cos ( —w)|): 
or le facteur de v 1 — e° est égal à 


2e[1+ecos(b—w)] — cos(* — w) [1 + e cos (P — v)]* 


=Í + ecos(9 — w]|2e — cos (# — w)[1 + ecos(?—x)]{. 


En vertu de cette valeur, le second terme devient 


m 5 [1 +e cos ($ — )] | [1 + e cos (P —5)] cos (# — w) (1 — y1 = &) -- 2ey1— e| 


e 


= — A [1+ e cos (P — w)] | [1 + e cos (£ — w)] cos ($ — w) — 4-2y1— etf. 
1+Vi—e 
|. On a donc 
moo EMO UE 25 2yr—e? . 
Rire [mac Mb e 


d'où, en mettant les valeurs (53) et (54), 


Tu s 
i dt 1+ V1— e 
(97) 
e cos( — w) 2y1—e? kX, 
x ————— + — — 0 
1+ V1— ei 1--ecos(o —v)| fu 


Les quatre formules (55), (56) et (57) étant supposées intégrées, on 
obtiendra les valeurs variables des éléments cz, e, a et e. Nous ferons re- 
marquer que, dans une orbite peu excentrique, la première équation (55) 
donnerait dz exprimé par une fonction ayant e en diviseur, et qui se préte- 
rait mal à l'intégratión. On aura de l'avantage à substituer à e et æ d'autres 
variables exemptes de ces inconvénients. 

Soient, en effet, 

t —ecoSo, 
(58) 


v = esins, 
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on aura 
dt Saa de AA e dw 
DT Nr dp" CRM 

58 br. 

( di do „æ de eds 
7 2— sin p + COST cd 


Substituant ici les valeurs (55), il viendra 
di B (N coss — M sins) + Asin, 


4 B (N sins + Mcoss») — Acoso. 


Pour effectuer plus facilement la transformation des coefficients de B, 
nous écrirons les valeurs de M et de N comme il suit : 


1 Erol 
M =2sin($ — 7) + 2" sin (2% — 2), 
1 3 
N = 2cos(9 — v) + " e cos (2 — 27) mz 
on en déduit immédiatement 


. I 
Ncoso — M sino = 2cos o + E ecos(24 — w) + 


NI 


ecosm, 


: PT epe A 
N sins + Mcoss = 2 sine + zesin (29 — v) + -e sinw. 


D | © 


Au moyen de ces valeurs et de celles de A et B, on obtient 


dio 2 COSP + Lecos(25 — v) 4 dictu TY k sin Xi 
di 2 DRE k fe 
uude 3 x 
(= [sine + Lesin(at —9) + 2 esins | Yir- kosb =. 
d 2 2 fe 


Comme on le voit, ces expressions sont affranchies de diviseurs qui pour- 


raient s'annuler. En les supposant intégrées, on obtiendra les valeurs de æ 

et e par les formules | 

e sins =v, 

(60) n 
ecosum = 6, 


d'oü z, sous la condition e 7» o. 
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Il n'y aura, dans l'emploi de ces formules, d'autre indétermination que 
celle de & quand e sera très petit; mais cela est dans la nature des choses et 
n'entraine aucun inconvénient. 

Détermination du plan de l'orbite. 


L'équation (11) donne, par voie d'intégration, 


(61) 9 — P + f'cos0dy, 
et les équations (18) sont 
do Zr 
(62) d. cosy, 
"o PENES UT 
(63) dtc gara sin e. 


Voici l'usage de ces formules : 

L'équation (61) fournirait une valeur approchée de o, en négligeant l'in- 
tégrale qui s'y trouve; alors l'équation (62) donnerait la valeur de 6, au 
moyen de quoi l'intégration de (63) ferait connaitre . On recommencerait 
les calculs en tenant compte de l'intégrale négligée la première fois, dans 
la valeur de 9. 

L'emploi de la formule (63) présente, dans le cas des faibles inclinaisons, 
le méme inconvénient qué celui de la première formule (55) dans le cas des 
faibles excentricités : il convient donc de recourir à de nouvelles variables 
auxiliaires. Posons, à cet effet, 


( F= sint cosy, 


(64) 


nous aurons 


| v' — sinô sin; 


dt d ; ; 
H” cosy cos 0 j,; — sin y sin0dy, 
dv» 


: d : 
um sin Ý cos0 VT + eos y sin0 dy; 


d'où l'on tire, en mettant les valeurs (62) et (63), et ayant recours aux 
relations (19) et (20), 


^ MR ze Ze Ar 
"e = — (cos0 cosy cos + sin y sing) — = — b FE 
(65) : à 3 
re t j vð loyer ILE 
m (cos0 sin cose — cost sino) xt b F" 
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où b et b sont les cosinus des angles que l'axe mobile des y, fait respecti- 
vement avec les axes fixes des y et des x. Ces équations étant censées inté- 


grées, on calculera Ø et Ý au moyen des formules (64). 


Il nous reste à changer, dans les formules (61) à (65), les signes de Ý et 5, 
de manière que ces quantités se rapportent aux usages des astronomes. 


Soient donc Q la longitude du nœud ascendant comptée, de l'axe fixe des x, 


vers les y, I l'inclinaison comprise entre o? et 180? : on aura ces autres 
formules 


p=% — f cosI d2, 


a + LT eos 
(66) ( dis k P 
m" BA. Ar Pa 
Sın di EE "u^ 9, 
! 1 ; : 
| = = — (cosI cos2 cose — sin 9 sino) ar e pes ; 
67 
(67) ES tat Q Ó si AS C. Zr 
FERÐ (cosI sin 9 cose + cos9 sinc) = = 
(68) sin] sin9 — + v’, 
sin Į cosg = — £', 


d'où Q et I; condition : sinl > o. 

Les derniéres de ces formules donneraient un résultat mal déterminé 
pour l'inclinaison I, si celle-ci était voisine de 90?; mais alors on peut déter- 
miner directement I par la deuxième équation (66) et tirer Q de la troi- 
sième de ces équations. Dans ce méme cas, l'emploi des formules (68) 
offrirait une vérification quant à la valeur de Q. 


Autres formes des expressions différentielles des éléments a, e, ©. 


Cette recherche n'a pas été faite par Wronski; elle a été suggérée par la 
simple inspection des résultats de la présente méthode. Suivant l'équa- 
tion (14), la variation de la constante des aires Æ est uniquement due à la 
composante Y, des forces perturbatrices. Nous nous proposons de faire voir 
que la variation du demi-grand axe a est due principalement à la compo- 
sante tangentielle des mémes forces, et que dans les expressions de e et z, 


5 
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la considération de cette composante, jointe à la composante normale dans 
le plan de l'orbite, permet de donner aux résultats une forme plus simple. 
Par le point M (fig. 2) de l'orbite, menons la tangente MT dans le sens 


Fig. 2. 


du mouvement de l'astre, et la normale MN dans le plan de l'orbite et 
dirigée dans le sens de sa concavité; l'axe des x, étant toujours dans le 
prolongement du rayon vecteur r, et l'axe des y, perpendiculaire à OM et 
faisant avec MT un angle aigu; nous désignerons par y l'angle X, MT. 
Désignant par db et ds l'angle décrit par le rayon vecteur et l'arc élé- 
mentaire de l'orbite parcouru pendant le temps dt, on aura entre ces diffé- 
rentielles, jointes à celle du rayon vecteur r, et l'angle y, les relations 


cosy ds — dr, 
sin y ds — rd, 
ou 


dr 
| PCOSY = a 


69 ( 
(69) : r db 
Ý SIN y = B : 


Projetons les trois composantes X,, Y, et Z, sur les directions MT, 
MN : Z, étant perpendiculaire à la fois à ces deux directions, si nous dési- 
gnons respectivement par T et N les composantes suivant MT et MN, nous 
aurons i> 

$ | T — -- X,cos1-- Y,siny, 
V) | N = — X, siny + Y, cosy; 
d'oü, inversement, 

( X,—Tocosq — Nsiny, 
(71) 


| Y, — Tsinq + N cos. 
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Or on a, en vertu de (27) et (45), 


dr fy 
scr Avi doo 
MT de k.esin( ©), 
2 
R rda- k, 
A CH à 


mettant ces valeurs dans les équations (69), nous en déduirons, sui- 
vant (71), 


| As dp: esin (t — w)T — fn]; 
(73) | 


Ny 2 E AE T pea sin (P — w) N]. 
La valeur dev sera d'ailleurs fournie par l'équation des forces vives (30), ou 


Esc MNA, 
e -^C i) 


relation qui, en y mettant la valeur de r en (6 — z), devient y = : E 
LI r^ 
{4 e eeg (P-co } 


Lum fo. [1 + 2e cos (P — s) + e]. 


Hn CUNSEWES Ote D : 


d'où, en vertu de la première équation (27), 
dE ——————— — — 
(74) ei accostu) e. 
t 


La substitution des valeurs (73) dans les expressions différentielles des 
éléments a, e et w étant effectuée, il nous restera à réduire les résultats à 
leur forme la plus simple. 

1° Demi-grand axe a. — Au moyen des valeurs (73), l'expression (56) 
devient d'abord 


1 da 2k 
A e(1— e) jiu + ecos(8 —w)] | 7— T2 < sin (b — w) N] 
ae esin (t — 5) Ë sin($ — w) T -75 N ng |: 
or le coefficient de T dans les | | est, en y mettant la valeur de ren (9 — x), 
hrec e 9I T esin? (P — Ð) _ : ks* "a(i1—e!) , 
k dut Jau = [re 2ecos (b— x) + ee DA k : 
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Le coefficient de N devient, en opérant de méme, 
I + e cos ( P — = 
A S 


I + eCOS(P — w) 
E 


e sin (b — v K 


on a donc cette expression trés simple 


(79) 


On voit ainsi que la variation du demi-grand axe est due à la composante 


| 


tangentielle T des forces perturbatrices 
2? Variations de l'excentricité e et de la longuude z du périhélie — Met- 


tons dans les équations (55) les valeurs (54), nous aurons plus explici- 


Lement 
e da á $ 1 
T = [2 + ecos(9 — w)] sin ($ — v) T OOE E A 


Ub. 


de 
= = [e + 2c00s(P — w) + e cos? ($ — w)] LL EE op LER n 


relations qui deviennent, en vertu de (73) 


Lt [2 + ecos(+ — «)] sine —m)| id | 
— cos(9 — a) | esin 0) T— 77. NT. 
T= s fle 2008 (8 — w) + ecost(& — æ)] | T+ 47 esin(#— w) N] 
+sin(®— w) [esin(#— 0) T — EN] 


Considérons la première de ces formules : le coefficient de T dans 
les | { se réduit visiblement à 
2sin ($ — 5); 
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le coefficient de N devient successivement 


fur. | k? 
er e sin? (P — c) a aco aea TT eor (e i), 


[2+ecos(b — »)] 


Mord t ipee ek (om) 


= e sin? (Ð — v) 
|. 2esin? (P — w) +e? cos (P — v) sin? (£ — 7) + cos (£ — z) + 2 e cos? (& — w) + e? cos? (P — v) 
E I + e cos (P — w) 


__2e+ cos($ — m) + e*cos( 6 — v) 
ae I + ecos(P? — v) 
| 2e[1+ ecos($ — w)] + (1— e?) cos (P — v) 


r ; 
= 2€ + — cos ($ — v); 
I+ ecos( — v) Ta E 


on a donc 


(76) LOREM 


dt asino — s) T+ [e 2 eos — m | N]. 


£ 
e 
Considérons actuellement l'expression de = nous aurons, pour le coef- 
ficient de T dansles| |, 
2 [e + cos (P — w)], 
et le coefficient de N deviendra 


e + 2c08( — m) + ecos? (P — m) 


sin (P — v 
( di 1+ecos(?— ð) 


— [1+ ecos(? —w)]} 


sin (b — 5) R 

= —— p— 2 2 E. a, i + bass 

Escola) t w) + e? cos? (b — w) —1—2ecos (P — w) — e cos? (P — »)] 
sin (P — 


I + e cos(P — v) 


ou, simplement, Ë 
— ; sin ($ — v); 

il s'ensuit 

(77) de (spei cos(® —m)] T— - sin(* NI 

77 di sl + CO ) — 7 sin( — a) {° 


Les résultats (75), (76) et (77) s'accordent avec les formules (122) de 
. ^ : À ; —Ó——— M — 
Wronski, quand on a égard à la différence des notations. - 
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Passons à la variation de la longitude moyenne de l'époque e. En vertu 
des valeurs (73), la formule (57) nous donnera 


Aot EEE e t o) [T A osito v] 
dt vi 1-- V1—e? 


= E ED | [esin m T- 77 Nji: 


IH Vie 1+ecos(P—w) hr. qd 


Le coefficient de T dans les | } se transforme comme il suit 


2 r : r I 
fun =) | ayi— e 1+Vi—e 


Le coefficient de N est 


1 2-recos(Ó —w) , | —— 
— | —— — esin’ (6 — w) + ecos (6 — v) [1 + e cos (£ — w)]! + 2y1— e? 
mer E ( ) ( ji ( 1] V 

|n (b — 5) + e cos (P — w) sin? ($ — x) | 
e cos (P — 2 e cos? (P — c) + e? cos? (P — w —— 4 
Ra E i $m) ( ) Tar 
I--V1— e? I--ecos(e — m) 
e 2e + cos( — m e? cos (b — m 
R l Da ( Len A 
1+ eCOS(® — ©) 
e e FROM SRE eec 
1+ V1— e? I +e cos (P — ð) V 


on a donc finalement 


` 


Fri r 
= = —1—2esin(d -—m)|[-————— — F 
dt "| ; j ayı— e? "Ne WIS 


1-5 e cos(bt — 5) +e \| 
+2 y1 e+ —— c= |e ——————— IN.. 
CeVicAN it ecos(b—v) 


\ 
Il resterait à transformer les expressions (59); mais il sera bien plus 


simple d'appliquer les formules (58 bis) aux valeurs transformées de e = 


x 
et - di? que nous écrirons d'abord comme il suit : 


de 1 x Pa Tem \ 
R = s (207 +2cos(#—w)T— ^ sin(# —a)N J 
edo 1 ia r ) 
T -s( Ch cp R MO NOMEN 
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On en déduit immédiatement, suivant les formules (58 bis), 


- NL ; ; Gg" : À 
-> — = | 2e cosw T + 2 cos T — - sinbN — 2esinwN | 
dt v a 
[449] I 4 H r * 
— — —-[2esinwo T + » sine? T + - cose N + 2ecosoN |); 
dt e a 
ou 
dt 1 P Á 
— = - | 2(cos& + ecoss) T — | - sing + 2esins | N |, 
i di. ...n ; a , 
(79) < $ 
| do I : ? r EPIS 
Ear 2(sin* + esinw)T + a CSP + 2e cosw N |: 


Ces formules sont plus simples que les formules (59). 


EXPOSÉ DE LA MÉTHODE DE WRONSKI ('). 


L'auteur de la Réforme absolue du savoir humain ne s'est pas borné à 
substituer, dans les problèmes de Mécanique céleste, les composantes des 
forces perturbatrices suivant les axes mobiles, aux composantes suivant les 
axes fixes; il en dispose autrement que ne l'ont fait les géomètres de son 
époque. | 

Les conditions qu'il s'impose produisent une simplification notable des 
expressions différentielles des éléments variables. Nous n'essayerons pas 
d'exposer les considérations toutes philosophiques sur lesquelles Wronski 
s'appuie pour établir ses équations; nous nous proposons d'arriver au 
méme résultat, en nous appuyant sur les principes universellement acceptés. 

L'idée qui a guidé Wronski est tout entière dans cette considération que 
l'attraction solaire et la composante X,, qui s'ajoutent algébriquement, doi- 
vent rester associées dans tous les calculs, comme elles le sont dans notre 


équation différentielle (16), sous la forme — = (fu — X,r*); nous avons 


2 


trouvé, par exemple, — pour expression du demi-paramétre (éq. 27). En 
pte, fü P P / 
k 


nous placant au point de vue de notre auteur, nous serons conduits à ajou- 


———————————— 


(1) Cette méthode constitue ce que Wronski appelle la Science de l’ordre. 
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ter à fu. la quantité — X,7?, ou à poser ya k k 
32 A 1⁄ 
k? kd 
(3a) P fe— X fr l 
^ 

Wronski considère en outre une quantité # ainsi définie Ws í 7 f 1 

gk k 
(81 ) W= C VT Se k j 


La premiere p est une ligne, la seconde, w, une vitesse, ainsi qu'il est 
facile de s'en assurer. Ces deux quantités p et w, qui sont fort peu variables, 
remplaceront le demi-grand axe et la constante des aires. Laissant de cóté, 
pour l'instant, la nature de ces quantités et les motifs de leur introduction 
dans les calculs, nous pourrons nous en servir, comme de constantes quel- 
conques devenues variables, pour arriver à l'intégration de l'équation (16). 


La combinaison des équations (80) et (81) nous donne 
pw = Kk, 


82 2 
va fu — Xi? = 2 cg. 


La constante des aires Á se trouve ainsi étre égale au produit des constantes 


p etw. 
Si nous mettons les valeurs précédentes dans l'équation (16), nous 
aurons 
ÆT cR vá Kk? 
det ciousm oo ge 


En multipliant tout par dr et intégrant ou indiquant l'intégration, il 


I dr* KS k? po “ul 
LE ns C kie Pd 2 Mov V — 98 0d unes 2 (2 Er 
2 Hp pr? T AE clar + fpo d r 3L M P 


En intégrant par parties, on aura 


I dr? w? 1 pw? "y TANT 
-LESCIET(iaeeifundge 


2 dt? 7 2 r* 


viendra 


soit actuellement 


(83) u=: f: d.p Í 5 d.p'w?, 
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notre équation deviendra 


dr* i apwe? l p'w! 


— — — J 


Mi F p? 
expression oü l'intégrale U est une quantité trés peu variable. 
Pour former l'équation différentielle de l'orbite, nous déduirons des 
relations (15) et (82) 


dE. pie 
dem p 


7 


En divisant membre à membre l'équation précédente par celle-ci, il 


viendra 
c dr? — ar | PES ug r? : 
rd p ` up 77 
soit, pour l'instant, 
(83 bis) p, d'oü Noms, 
> ^ p x 


s dr D AA g 
nous aurons, en substituant la valeur de . et celle de £, dans l'équation 


précédente, 


ou, en multipliant tout par p? d?, 


j dp V? U^ U 
(84) (4 P) = (ae — e — is) der 1-5 (=? dé. 


w? 
Or, si Pon fait abstraction de la variation des quantités p, U et w, et que 


l'on pose 
U 


(85) I > 


on trouve facilement que l'intégrale de l'équation (84) est 
(86) =1 + e cos( — m), 


z désignant d'ailleurs la constante d'intégration. 
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En ayant égard à la valeur (83 bis) de P, on obtient, pour équation de 


l'orbite, : 
^ -— T AR CMS LI 
(87) der I + ecos(* — w) 


Il s'agit actuellement de disposer de la variation des constantes, de ma- 
nière à satisfaire à l'équation différentielle (84), qui est du premier ordre 
seulement. Nous aurons, en différentiant (86) et y faisant tout varier, 

(87 bis) do = cos (Ð — w) de — esin (P — v) d ($ — w): 

mettons cette valeur dans (84) et observons que son second membre se ré- 
duit, en vertu de (85) et (86), à e?sin?(6 — z)d$*; cela nous permettra 
d'extraire les racines des deux membres et nous donnera 


cos(P — »)de —? T — esin (P — w) (d^ — dw) = — esin(* — v) d^. 


Nous n'écrivons pas le double signe, attendu que, notre équation devant 
encore étre satisfaite dans le cas oü les forces perturbatrices seraient nulles, 


le signe — remplit seul cette condition. 
Mettant à la place de p sa valeur P, et supprimant les termes en d6 qui 
se détruisent, il restera 
(88) cos(* — w) de — Ë + - sin(P — wv) edo — o. 
A cette équation de condition va se joindre celle que nous déduirons 
de (85), en l'écrivant ainsi 


U 
I — e? m E 
d'où, en différentiant, 
aede. dU. adw | : dU 2 div 
1—e . U Ww. (te We 
ou 
d 1 dU 
cde jio qe ERN 
w 2 w 


Or, nous tirons de (83), 
2 2 
dU = x (pd.w? + w* dp) — a "pd.wp, 


1 dU 1 dw dw 
EU V = (pZ dp) 5 (rZ + dp); 
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Mettant cette valeur dans celle de ede, on aura 


N dp 2p ?N dw 
|æ AERE: em 
/ \ E 7 wW 


r P 


ede P (£— 


Or on a, en vertu de l'équation de l'orbite, 


£ 


R = ecos(9 — w); 


d'un autre cóté, le coefficient du second terme peut s'écrire 


(£ — iJ —e=— esin (9—0); 


p 


substituant ces valeurs et divisant tout par e, il vient finalement 


l 
( 89) de =£ cos($ — w) P — esin? (b — a) Z. 


Au moyen de cette valeur, l'équation (88) va nous fournir celle de eds: 
elle devient 
div 


lp 1 
P cost (b — w) T — esin'(& — w) cos(b — w) — — £ r FR En) edu -3o; 


d'où, en réduisant et divisant ensuite par sin(® — >), 


4 1 p ap dw 
(90) e dw = Sin (Ð — w) ») |£ = + ecos(9 — v) — — 
Cette expression peut étre mise sous une forme plus simple : on a, en 
effet, 
dw d.pw dp 


(91) abe oe bey. 


la substitution de cette valeur conduit à la formule 


(92) edə = sin (P -— v) | Æ +ecos(# — d.pw À 
; P pw 

On aurait pareillement 

(92 bs) de = [cos ($ d cepe > — esin? (ee) mm. 
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Nous obtenons ainsi (89) et (90) ou (92), les variations différentielles de 
l'excentricité e et de la longitude z du périhélie, en fonction des variations 
des constantes p et w, dont les valeurs explicites deviennent, en substi- 
tuant l'expression (14) de £, dans (80) et (81), 

K (f Y,rdt dt)? 

| fre Ju — Due WU 


USA iue c Ud 
T f Yirdt 


(93) 


La première de ces expressions fait connaitre le demi-paramètre. Nous 
avons donc trois des quatre expressions différentielles qui sont nécessaires 
à la détermination du mouvement dans le plan de l'orbite. Quant à la qua- 
trième, Wronski n'a pas à s'en préoccuper, attendu qu'il détermine la longi- 
tude ® par l'intégration directe de l'équation (15), qui devient actuellement, 
en vertu de (82), 
de = ET qt, 


r 


(94) 


On fait ainsi disparaitre les difficultés auxquelles donne lieu la variation de 


la longitude moyenne de l'époque. 


Nous avons obtenu la valeur de r (87); il nous reste à obtenir sa dérivée 
en fonction des nouvelles constantes et à calculer la vitesse e dans ce 
systéme. 

Nous tirerons la valeur de 7 ; des équations (83 bis) et (87 bis), ce qui 
nous donnera 

| pa . " Á 
— = — — -[cos(* — v) de — esin(* — w) db + e sin ($£ — v) dw] 
r P P 
ou 


e enitn Cou) di - ; | c08( — m de — P P + sin(s— s) eds]. 
d P P r p 


expression dont le deuxième terme est nul en vertu de (88); on a donc 


simplement 
dr = esin(® — w) d, 
P 


expression identique avec celle que fournirait l'équation (87), différentiée 
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' , HEP, dr K 
seulement par rapport à ®; en d'autres termes, la dérivée jg conserve [a 


méme forme, dans lé mouvement troublé, que dans le mouvement ellip- 
tique. 

Mettons pour 7? d6 sa valeur £dt, et remplacons 7 par w, suivant (82), 
nous aurons 


dr q 
(95) Ji = wesin (P — v), 


et l'équation (91) nous donnera 


- rd% p 
6 E nsu 
; z dr de UAT i 
La forme des expressions de 3 et R est ainsi la méme que dans le 


mouvement elliptique : ces valeurs ne different de celles de l'ancien svs- 
tème qu'en ce que, dans les expressions (80) et (81) de p et de w, la con- 
r a 2 IET ERR à ia 
stante de fp se trouve accompagnée de la quantité — X,7?, qui n'existe pas 
dans les formules de l'ancien système. 
En élevant au carré les expressions précédentes et ajoutant, on a 
/ p? 
g^ — wt (Z + esinta), 


r? 


expression qui se transforme en 
(97) 9*— w [1+ 2ecos(* — v) 4- e]. 

Considérons les deux valeurs e, et e, de e, qui répondent au périhélie 
et à l'aphélie : cos(® — æ) étant, dans le premier cas, égal à +1 et 
dans le second égal à — 1, l'équation précédente donne e, = œ (1 + €), 
v= w(1— e); on en déduit 


I 
$9 — » (oi + £3), 
relation qui précise la signification de la constante ø. 


Appliquons actuellement les équations (58 bis) au calcul des dérivées des 
auxiliaires Ë et v, dont les différentielles peuvent être utilement substituées 
à celles de e et æ: nous aurons, en faisant usage des expressions (89) 


WWwW.rcin.org.pl 


B.46 YVON VILLARCEAU. 


et (90). 
eos i d 
| PN SNC en esint sin (t — a) —; 
x p w 
(98) dw 


à RE: | À 
| dv =Ë sine © -+ecosPsin(P— w) — > 
| g P w 


expressions beaucoup plus simples que celles auxquelles donne lieu l'em- 
ploi de l'ancienne méthode. Leur intégration étant effectuée, on aura, pour 


déterminer e, z, les relations (60), 


(99) condition : e > o 


Comme on le voit, la détermination des éléments variables, par l'emploi des 
constantes p et w, est d'une simplicité trés remarquable et semble devoir 
motiver la préférence en faveur de la nouvelle méthode. 

La détermination des constantes qui fixent la position du plan de l'orbite 
se fera suivant les formules (66) à (68). 


Pour compléter le système des équations de Wronski, nous devons 


de 


ajouter une nouvelle détermination dont voici l'origine. La relation r — T 


(p. 15) suppose nécessairement l'angle ® et, par suite, z, comptés d'une 
méme ligne fixe dans le plan de l'orbite : si donc on désigne par x la lon- 
gitude de cette ligne fixe, comptée suivant les usages des astronomes, en 
sorte que l'on ait, pour la distance au nœud, cette expression 


$-—y--9-—8&2, 


on en déduira, suivant (66), 


y — Q2 = ọ — b — — fcosldo 


ou 
d(y — 9) = — cosId9, 
puis 
- d y — Q 4 A 
(100) sin] a = — z4 sino coslI. 


Cette expression coïncide avec la formule (125) de Wronski, moyennant 
une correction à appliquer à cette dernière, comme aux équations (70) de 


ses Prolégomènes : cette correction consiste à multiplier les éléments diffé- .- 
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rentiels par le facteur 


Vi+ 2e cos o + e? 
I + € COS o 


facteur qui diffère peu de l'unité, quand l'excentricité est trés petite. 


ee d(y — 9 k : : : : 
La valeur précédente de en devient infinie dans le cas d'une incli- 


naison nulle et ne se préte pas à l'intégration dans cette circonstance; il en 


fc» Sa MD Á : EXT Å 
est autrement de la dérivée Ý : en effet, si l'on ajoute la troisieme équa- 


tion (66) avec (100), on a 


sinl dx = (1— cosT) ^^ sinc; 


dt 
d'oü 
(101) x = = tang +I sins. 


On a dü remarquer que le demi-grand axe ne figure pas dans les équa- 
tions différentielles, propres à la nouvelle méthode; il n'est pas sans intérét 
d'en effectuer la détermination, ne füt-ce que pour montrer la grande diffé- 
rence de son expression avec la formule (56). 

Les constantes du mouvement elliptique sont liées par la relation 


SES wt 
PF T. 
d'où 
d 
da _ dp 3e de, 
a p I 
et, en vertu de (89), 
2 , 
(102) ee | ee i E cos(+ — w) | Aya ; sin? ($ — s) Ld 
a V= gg P Fee” w 


Au moyen de l'expression (92 bis) de de, on aurait cette autre valeur 


da 2e[cos(P— v be n l. piv 
da .yae[cos(b —m)--e] | || dn... 28 as uy d- Pr", 
a ( 1— e? Cp 1 — e? pw 
ou 
da I+ 2e COS(P — v 2 d . piv 
(103) — Æ Vi scribe puits Jam a" VERAM 5 sin? (^ Ms w) d pw 
a Bið P 1 — e? pw 
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expression qui, en vertu de la relation (97), peut encore s'écrire 


n4 2 

et dp 2e ^ d. pw 
ARA SP a ER sin? (P — D) n . 
e*(r—e?) p 1— e pw 


I 
(104) z da= 


Si l'on utilise cette formule pour le calcul de a, on pourra lui laisser 
la forme (104), qui est composée, comme celle (92) de edz, au moyen des 
différentielles de p et de pw. Mais, pour effectuer la comparaison que nous 
avons en vue, il importe d'exprimer explicitement la valeur de da, au 
moyen des composantes des forces perturbatrices. 

La premiere équation (93) nous donne, en vertu de (14), 


dp -: dk | diXr*y 


p ET Ju 
on a d'ailleurs, suivant (82), 
d.pw _ dk. 
pw cs 


la substitution de ces valeurs dans (104) nous donne, en ayant égard à (97), 


dk 1+2ecos(P—m)+e? d(X,r?) 
> Es am — — , 
[1+ e cos( — w)] rM = fara 


1 
~= da = ; 
a 1 — e? 


ou, en mettant pour d£ sa valeur Y,rdt, et transformant le premier terme, 
1 


Yi r1--2ecos(P$ —w)--e* d(X,r?) 


I et un i : 
z da = (4 ecos( $)] i dt + TEE T posso 
puis, en vertu de (82), 
I FE 2k Yidi 
[09 ë 
pe 3 12-2ecos(b—9)--e* d(X,r*). 
1 — e? fu—X;r? 


Sous cette forme, on voit que le premier terme de cette expression 
diffère, du premier terme de l'équation (56), par des quantités du second 
ordre relativement aux forces perturbatrices; quant au deuxieme terme, si 
l'on développe la différentielle, on aura (95) 


d(X,r?) —2X,rdr + r*dX, =2 X, rwesin (t — w)dt + r!dX,; 
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7 k 4 " . 
en mettant ici pour w sa valeur = le deuxième terme de (105) deviendra 


CURE 
A sin (£ — v) 


1--2ecos(P —5)--e?  X,dt. 
I + € cos($ — ð) fv»—X,r 
Moka 2ecos(b—m)+e r dX, 


1— e? fu—Xir 


Le premier de ces termes diffère du terme en X,, dans (56), de quantités 
qui sont du second ordre par rapport à l'excentricité et aux forces pertur- 
batrices. Quant au terme en dX,, il n'en existe pas de cette nature 
dans (56). 

Telles sont les différences que présente la quantité da dans les deux 


méthodes. 


Bien que Wronski n'emploie pas la variation de e, il est cependant inté- 
ressant d'en former l'expression dans le système actuel. 

Pour cela, nous déterminerons d'abord le moyen mouvement. L'équa- 
tion (95) nous donne 


dr? 
dt? 


1 


2 
= w e? sin? (P — v) = w? [e — e? cos? (6 — w)] = w? |e -= (: — ) Ë 
d'ou, en multipliant par 7°, 


dr? ; 
r? = = w? (er — r+ 2pr — p) — w*(r— ey) [ee — (r — ay*] 


dt? 
dr* = 2 
ph — —w'gt(r—e* Mei q 
de ( | NM: )) | 


Nous intégrerons cette équation en négligeant la variation des constantes 


ou 


et posant 
+ 
(106) — — 1 = e Cosu, 
a 
d'oü 
| dr ed du 
=- = Æ — e Tr; 
a dt dt? 


substituant ces valeurs dans l'équation précédente et supprimant les fac- 


7 
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teurs communs, on aura 


du? ` «æ?(1—e?) 
2 EE . 
(1— e cosu) ET OR TT 


d'où, en extrayant les racines et posant 


VE 
(106 bis) Ive 


a 


(1— ecosu) du — ndt: 


nous avons choisi le signe qui donne à du le signe de dl. 
On a, en intégrant et en désignant la constante par £ — z, 


(107) fndt -- « — = u — esinu. 


Il était nécessaire d'effectuer l'intégration pour obtenir la valeur der 
dans le nouveau systeme ou, du moins, sa relation avec le demi-grand 


axe a. De la relation (106 bis), nous déduisons 


w? 
a? n? — w*(r— e*) — d 
a 
ou, en vertu de la seconde équation (82), 
(108) a? n? — f y. — X,r?, 


relation à laquelle on pouvait s'attendre. 
Pour avoir égard actuellement à la variation des constantes, différentions 


complètement l'équation (107), nous aurons 


de. do ` í du E 
(109) atg Tg 07 e0088)-,, — RUES 


d'un autre côté, l'équation (106), mise sous la forme 


* 
— —1—eCO0Su, 
` Ø 


nous donnera 
dr da _ esinudu — cosude 


— — — LLL ma 3 


r a I.— e COS U 
ou 
r {1 dr 1 da ne TIR du dux de 
Gave dt uuu dt dt 
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Ei A du : $ 
Eliminons a entre cette équation et (109), nous aurons 


dr D= de 
E c Iq (cosu — e) F: 


Considérons le terme en Æ - 


z? Ce terme, en y mettant la valeur (95), est 


7 A we P . 
ttes we sin (D — 5) = — t /1— sinu — — nesinu; 
ce qui détruit le terme en 2, comme cela doit étre : il reste 
e T QE sinu — sE + (cos u — e) m 
dé. da ws dt dt? 


d'où, en vertu des relations (40), 


(eds — edw) sin (b — v) =— y 1 — e? = = + V1 — ecos (P — w) de. 


Au moyen des valeurs (90), (102) et (89), on en déduit 


= 


ed: sin (b — z) = sin? (b — v) [E2 E + e cos (P — D) "s 


qp. um ; sin? ($ CARTES e zs 
P Læt 


== 2e 
Agree pq 865 P angi. "2 
all 1— er 


d 5 d 
+ Vie cos m|^ cos (D — m) £ —esin (8 — $2 |. 


Occupons-nous d'abord des termes en «1 — e*; nous aurons, pour ceux 


en dp, 
teur ad 2(q — As A i dp. 
Vin et € Ë cos? (4 7) 2ecos (b — o) | P p 


Les deux termes négatifs renfermés entre crochets étant réduits au même 


dénominateur sont 


I — e? + 2ecos($ — w) + 2e? cos (b — a) — e sin (P— v) — [1+ e cos(t —»)J*. 


2 


di 1 + e cos (Ð — D) I + e cos (P — s) 
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le terme dont il s'agit peut donc s'écrire, en mettant pour P sa valeur 


Vie 


za COS (£ — ø) [124- e cos (4 — m) ]22- esin? (9 — 5) — [1+ e cos( — a)]? dp 
1 + e cos ($ — D) P 
K | esin? (b — w) — [1 + ecos (b — æ)]? sin? (£ — 9)j ap 


—wVi1—et E sin? (D — v) [1+ ecos ($ — 5) — e| a 
P P 
Quant aux termes en y1 — e* c ils sont réductibles à 


Vi— e — e? sin? (P — v) JEP — cos (^ -|F dw 


On a donc, en supprimant le facteur commun sin(® — z), 


ed: = sin (£ — v) |? a + e cos (P — x) = 
rp w 
— Vie sine — a) [+ ecos(# m e| 7 
: dw. 
+yi— ee sin — ea neZ — cos( —9 |: 


d 3 + ; 
or les termes en Z se combinent de la manière suivante : 


eta d D eaae el oT Frau dp 
sin ( IC Vi— es + eyi E: 
=sin(# —a) (7 


dp 
2,/1-— ei ` 
Ps te | ne 


dw ; : 4 
Les termes en v ose réduisent à 


esin (9 — z) | cos(# — 9) — cos 8 — 9) Sea de aie r | Æ dad 


w 
= e sin ($ — x) e: 


dw 
--9evi—e? 
ues E 


w 


Mettant ces valeurs dans l'équation précédente et supprimant le facteur 
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commun e, il viendra 
Ë ! —= a : 
d:—  esin(b—5x) pico eee as 5m 


ENADUA i 
ES I 4- V1 e P 


n 


t esine — o) | ai ar + 


e dw 
— COS (P — 5) E 
P I--VI—e* 


En introduisant pw à la place de w, on aurait cette autre expression 


pa sine — w) | horn cars 
Da e PiP 
t 9 
| estes): pego d. £ cos (P — m) p 
P 1+V1—e pw 


Nous ferons remarquer que les valeurs des différentielles de, edz et de, 


TIP d. pw n "De j 
exprimées au moyen de —P*., sont au moins aussi simples que celles expri- 
pø 


d.pw _ dk _ Yırdt 
DIS Ski: ^ 
l'avantage que Wronski a pu trouver à substituer la constante # à la con- 


å dw i Ë 
mées au moyen de = Or, comme » on ne voit pas 


stante k. Notons, toutefois, que, dans le Tableau final de ses formules, il 
n'emploie la constante w isolée, que pour le calcul de e. 


Nous terminerons cet exposé de la méthode de Wronski en rapportant 
ici, sans les démontrer, les deux formules suivantes : 


pe : 
I + e c0s(P — m) 


(113) == f(1— sin? ($ — w) Er d.piw 


1--ecos(P—o)/ p 1 + e cos (P — w) pw 


(112) koe 


oü la lettre e désigne la base des logarithmes népériens, tandis que p, 
et æ, sont des constantes égales aux valeurs que prennent les éléments 
variables p et æ à une époque arbitraire et donnée : 

La constante de l'intégrale € est censée déterminée de manière que Z 
sannule à l'époque donnée. 
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Wronski profite de ce que les éléments du calcul des intégrales qui figu- 
rent dans l'expression de X se trouvent tout préparés, en vue d'autres dé- 
terminations, pour obtenir une expression relativement simple du rayon 
vecteur r, en substituant l'intégrale E dans la formule (112). 

[l est facile, étant donnée la forme (112) du rayon vecteur, de rechercher 
la valeur de l'exposant €, qui est nécessaire pour tenir compte de la varia- 

ion des éléments e, w el p. Nous avons fait ce calcul, et nous avons trouvé 

une valeur de € dont l'expression (113) contient seulement les termes 
principaux. On nous dispensera de reproduire ici ce calcul, si l'on observe 
que la valeur complète de £ sera nécessairement plus compliquée que l'ex- 
pression (113), et, par suite, dépourvue des avantages que l'on pourrait 
songer à utiliser dans le cas d'une forme plus simple. 

Au reste, Wronski parait avoir reconnu lui-méme les inconvénients de 
l'emploi des formules précédentes; car, aprés avoir reproduit, sous la 
marque (197), la désignation des intégrales qui constituent, suivant lui, la 
science de l'ordre, il s'exprime ainsi : 

« Mais il est inutile, pour la détermination de la quantité de E, d'opérer 
la dernière (178) de ces intégrations, parce que, comme nous l'avons déjà 
dit, la solution du problème-universel, dans laquelle entre cette quantité €, 
ne doit servir que pour des cas spéciaux, oü, d’après ce que nous avons 
indiqué sous les marques (74) et (75), cette intégration peut suffisamment 
être opérée par les nombres constants p. et y. Pour l'exacte solution du pro- 
bleme-universel, il faut toujours employer la formule (181) (expression ordi- 
naire du rayon vecteur), qui ne contient plus aucune intégration immédiate. 
Il ne reste donc qu'à opérer les sept premières de ces intégrations (197) 
(voir ci-dessous), pour achever complètement la théorie lunaire, et géné- 
ralement toutes les questions pratiques de la Mécanique céleste. » 


Résumé de la méthode de Wronski. 


Le principe fondamental de cette méthode consiste à combiner, avec l'at- 
traction exercée par le Soleil sur la masse m, la composante X, des forces 
perturbatrices suivant le rayon vecteur r. 

k désignant la constante des aires, p le demi-parametre et øw la demi- 
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somme des vitesses au périhélie et à l’aphélie, on a, entre ces quantités, les 
relations 


W= 


A k? TESK 
P= før? E uu. 


pw -— Kk, 


et l'on peut prendre arbitrairement l'une de ces trois quantités comme 
élément distinct de l'orbite. 

L'excentricité e, le demi-grand axe a et le demi-paramètre p conservent 
la relation 


p — a(1— et), 


tandis que le demi-grand axe et le moyen mouvement z sont liés par cette 


nouvelle relation 
ær sf —X,rh 


Wronski conserve l'excentricité et il substitue à la longitude du périhélie 
celle de l'aphélie; nous ferons usage de la longitude ø du périhélie. Quant . 
au quatrième élément nécessaire à la détermination du mouvement dans le 
plan del'orbite, notre auteur substitue à la longitude moyenne e de l'époque 
la longitude vraie à cette époque. 

Enfin, la situation du plan de l'orbite reste fixée par la longitude du 
nœud ascendant Q et l'inclinaison I. Å ces éléments s'ajoute la longitude 7 
de la ligne fixe dans le plan de l'orbite, d’où se comptent les angles ® et z, 
de telle sorte que la longitude du rayon vecteur et celle du périhélie, 
comptées suivant les usages des astronomes, sont respectivement 


Ate et x-rw(!). 
Ceci posé, voici les expressions des diverses constantes variables dont il 


est fait usage dans la nouvelle méthode : les formules sont classées sous la 
marque (197), p. clxviii de la Réforme absolue du savoir humain. 


(1) Entrainé par l'habitude, il a pu nous arriver d'appliquer aux angles 9 et z la dénomi- 
nation ordinaire de longitudes; mais le lecteur aura sans doute suppléé à ce que ces déno- 
minations peuvent avoir d'incorrect. 
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Tableau des formules de Wronski. 


Numéros du présent Numéros 
Mémoire. de Wronski. 
(Fi yet4Ba). .. LB) dices. pw=k= f Yirat, 
1 tp y. dp dw 
(89) A Ts e = P cos(e — a) P f e sin? (Ré mn) ss 
{ » 


UL z 
(66) et (99 bis). 1'*éq.(70)... 1l AL x cos(y, — 9 + D) dt, 


(1013... «s 9*6]. 20) s = 
+ > k pw 
Hekk 75 Eos me ou Li ia dt, 
£ T ey ePi 
(GA VT (Er). AN msn P) P cos(a — j 
pie. 
(66) et (99 bis). 2* éq. (70)... o= f sin (y — 9 + 9) —. 
9 j k cl 
Ecl Fo Iu) e: 
ye 1--ecos(P — 5) | p 


e cos (6 — 5) sin? in — 5%) d. d.pse 
I + e COS (P — w) pw 


(Lais ieu aðar CE DT sr 
| (5: 


Cette dernière intégrale sert à calculer le rayon vecteur suivant la for- 
mule 


poe 


(ria. oo 9 (180). ...... "= 1-recos(9 — m). 


où p, et ø, désignent respectivement les valeurs de p et æ correspon- 
dantes à & — o. 

Nous ajouterons à ces formules quelques-unes de celles que nous avons 
démontrées, et dont Wronski ne parait pas devoir se servir : 


(1) Se reporter à la citation du texte de Wronski, reproduit page 54. 
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Numéros du présent 

Mémoire. 

à ` dp Á d.pw 
(oa bts)... 1 , e= | [cos(® — s) + Le — | esin?($ — x) 


pw 


? 


| EL P cos P — Sesine sin —9) $, 

OL NE d i d 
| = (P sine P + fecosæsin(æ — o) Z, 

r p w 


ot les deux suivantes qui se déduisent aisément de (98), 


k d EH ; d.pw 
ë= f (cost + ecosm) P — [esine sin(# — n) E 
P pw 
(y 3 1: ^ MSRP , , 
pw 
| 2 (sing + e sins) Z + f ecosæsin(æ — n) í EE 
P pw 
e sinu =v, 
(HOVER im. | a 
€CoSo = 5; 
d'où æ et e. 
Log bre =y — 2+, 
I 4 À Zr 
| ë= — f (eost cose cose — sin9 sine) E dt, 
NOTES ets : Lr 
| =+ f (cosIsine cos e + cos sing) 7^ dt; 
ou 
; 1115-5 ab dap 
Ë =— [| costx+ 5) — acos cosg sin? = 1| S^ ar, 
TIU A i 
£e f [sino P) — 2 cose sin9 sin? TIE dt; 
{ sinI sinQ— + v', 
(OMR. j 
sin Icosg= — ý’; 


d'où Q et I. 
Dans la Note suivante, nous donnerons les expressions des composantes 


X,, Y,, Z, des forces perturbatrices. 


DELL ooo, A RR 


ð 
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CALCUL DES COMPOSANTES X ,, Y,, Z, DES FORCES PERTURBATRICES, 
SUIVANT LES AXES MOBILES; 


PAR M. YVON VILLARCEAU. 


Désignant par m’ la masse de l'une quelconque des planètes troublantes, 
æ', Y', z’ ses coordonnées rapportées aux axes fixes, z’ son rayon vecteur 
et s’ sa distance à la planète troublée m, on a d'abord 


(116) $'? — r* — 9rr! cos(r, r') 4- r'*. 


in ne considérant que la seule planète perturbatrice m, on a pour 
expressions des composantes, suivant les axes fixes, 


I / 
L — X c 
S 1 
x = fm "E 5) 


7 


Or, si l'on fait usage des notations (1), on aura, pour expressions des 
composantes suivant les axes mobiles, 


X,— aX -- a' Y +a", 
Y,— 0X + 0/ Y + b'Z, 
Z,-—cX-r-c'Y--c2; 


mettant ici les valeurs précédentes de X, Y, Z, il vient 


fm 
Kar 


fm 


[a(x'—æ)+a(y—y)+ a'(2 — 2)] — 5 


(ax'+ a'y'+ a" s'), 


Y,= I [bla — 2) + (y! — y) + ba — 2)] — ER (oa + d'y + ba), 


Z, = I [e(a! — 2) + e'(? — y) + e'(! — 2] — 05 


PE (cx'+c'y +c" a ; 
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expressions qui peuvent encore s'écrire 


Å F — 5 )(azt- a! y! -- a" s) — TE ars a! y -- a" 3), 
ID) V = fm ( 5 — 5 (zt b'y'+ b^ s!) — 2 (ba + b'y + 0" z), 
| Z; = 1m (> — zat c! y! 4- c" 3) — I (ez + c'y + c" a). 
Or, on: a 
VET conr. ow), asp cog (7, 5); m enods(1 ka): 


il s'ensuit, en ayant égard à (1), 


| a x' 4- a! y' 4- a^ 3 — r' cos(r', 24), 
(118) € bz! 4- D' y! + b z — r' eos (r', y), 


| Cz -- c yt 4- e a r'eos(r, z,). 


On aurait de méme 
ax 4- a' y + a'z —rceos(r, xi), 
bx b'y =+ b"z =r cos(r, yi), 


cæ+c'y+c'z=rcos(r, 3). 


Le premier des derniers cosinus est égal à l'unité, et les deux autres sont 
nuls. En vertu de ces relations, les formules (117) prennent cette autre 


forme 
i " fm 
X,—/ms P zi) r'cos(r', 21) — : ys 7? 
(119) Y ca = M cos (7, yi) 
c ` s” r EE 


I I 
d = fm 3 — =)" eos ( r', 31). 

Le calcul des cosinus qui restent ici peut se faire de deux manières diffé- 
rentes : 1° au moyen de la seule trigonométrie sphérique ; 2° en faisant 
usage des expressions (118). Nous nous bornerons à l'emploi des for- 


mules (118). 
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Avant tout, nous ferons remarquer que, l'axe des z, coincidant avec le 
rayon vecteur r, le cosinus qui figure dans l'équation (116) a pour expression 


(120) cosi r, Free e08(7, d). 


Calculons d'abord les coordonnées +’, y’, z en fonctions des éléments 
de l'orbite de m. 
Soient Q’ et T la longitude du nœud et l'inclinaison de l'orbite de m, 


Fig. 3. 


ss M 


7 la distance de r’ au noeud N'; on aura ( fig. 3) 


ato p sos ri. Ah Fr ORC S N z' — r' cos(r', Z) 
ou 
| æ'= r'(cosQ' cose' — sin 9' sine’ cosI'), 
(121) { y'= r'(sin 2’ cosy’ + cos’ sing’ cosT'), 


L2 irait sini. 
Changeons maintenant, dans les formules (5) de la page (64) du 2° Volume 


de la Mécanique de Poisson, 4 en — Q et Ø en — I, nous aurons 


a! — + cosI cos9 sinc + sin 9 cos c, 


7 


a" — -- sin I sin o, 


b — — cosI sin 9 cosy — cos? sin c, 
(122) ( b'—+ cosIcos9 cose — sin 9 sin o, 


b” =+ sin I cos, 


e —-- sinIsin 9, 
c' = — sin I cosg, 
e= E cosi. 
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Au moyen de ces valeurs, on obtient 


a x' = r'(— cosI sin 9 sinc + cos9 cos) (cosQ' cosc' — sin 9' sinc! cos’), 


a' y'= r'(-- cosI cosg sine + sin 9 cos?) (sin 9" cos ' + cos9' sing’ cost’). 
Effectuant les produits et ajoutant, il vient 


| + cosI sine cos«' sin(9' — 9) + cose! cose cos(9/— Q) i 
a x'-r a! y! =P : : Å i 16 
j — cosT' sin«' cos sin(2'— 9) + sin sin ọ' cos(9/—90) cosI cosl | 


Or, nous avons 


" I 
cosi — 1— 2 sin*- I, 
| deed t -— 2! I 
cos = 1 — a sin SE 


z I A I A I á I 
cosI cosl'— 1 — 2 sin? 1 — 3 sin - I 4 sin? z sin? =Y. 
2 


Substituant ces valeurs, on trouve aisément 


/ i 


| + Sin 9 cosy sin(9'— 9) — 2 sinc cosy’ sin (Q'— 2) sin? ; I 
— sin v cose sin(9/— 2) — 2 sinc sin g' cos (9/— 0) sin? = I 
ax'+ a! y'= r'{ + coso cose! cos(9' — 9) + 2 sing’ cosy sin (9/— Q) sin? Á L ) 
+ sin sing’ cos (2'— 2) — 2 sin ọ' sin y cos (Q'— 2) sin? ; V 


PINE 
i in 5» C 2:1 À 
| + 4 sino sin + cos(9'— 9) sin* - I sin =t 
ou 


* . * . * 9 1 
— sin (?'— 9) sin (Q'—Q) — 2 sing sin(s'-- 9' — 9o) sin? I 
à i Fg] 
ax'+ a! y'= r' | + cos (9'— 9) cos(0'— 9) — 2 sine'sin(g + Q — 0^) sin? 3 r 


| 
| 
/ 
` è I è 1 152 I r 
+ 4 sino sino cos (9/— 9) sin? = I sin? — 


Les termes de la | qui sont indépendants des inclinaisons se réduisent à 


: cos [e' + 9' — (o + 9)]. 
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On a d'ailleurs 
a" z' — r' sine sine' sin I sinf’. 


En ajoutant ces termes, pour former la valeur de cos(r', x,), suivant la 
premiere équation (118), nous en transformerons les éléments, d'apres la 


relation (99 bts), 


(123) p—y+b—Q ou o--Q-—y-4- 9, 


-6 


et nous poserons 


|= 


«= — sin (y + 9 -— Q) sin(y + Ð/— 9) sin? - I 


T 
4 > 


: ` » I 
— sin(y/4- 9?'— Q) sin (y + ? — Q) sin? = r 


N 


, Å E EE 
+ 2sin (y + P — Q) sin ( y/ + 9/— 90^) cos (Q'— Q) sin? X I sin*- à 


Los t ; à S 
^j sin(y + * — 9) sin(y^4- $' — 9^) sinI sinf’. 


De cette manière, nous aurons (118), (119) et (120), 


à dug 
(125) r'cos&(F,r) = TER oem | cesar Ð'— y — P) + 2 «| 


Il est facile de déduire de ces formules celles qui expriment la valeur de 
r'cos(r', y,); en effet, pour substituer l'axe mobile des y, à l'axe des æ,, il 
suffit d'augmenter la distance 9 au nœud, ou l'angle 9, de 90%. En consé- 
quence, si nous désignons par B ce que devient g, quand on augmente ® 


de 90?, nous aurons 


1 ; v 
> B = — eos(y + P — 9) sin(y'+ 9'— 2) sin? — I 
4 ) 
: x od 
— sin (y/4- ?'— Q^) cos(y + P — Q^) sin? ^ rU 
(126) ! 
; st es 74 
+ 2C08(y + P — 9) sin (y' + *?'— Q^) cos (9' — 9) sin? - I sin? = d 
I À á P 
m COS (y + P — 2) sin(y'4- ?/— 9") sin I sinl’, 
. I 
(127) r'eos(r*, y) =r | singe 8— y — e) 18]. 
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Passons au calcul de la troisième équation (118). Des formules (121) et 
(122), nous déduirons 


cx 4- c! y! -- c" z' — r'[+ sin] sing (cos9' cose'— sin9'since'coslI) 
— sin I cosQ ( sin 9' cos«'4- cos9'sin«'cosI') + sinI' cosIsins'], 


ou, suivant la troisième équation (118), 


r' cos(r', 21) —r'[— sin I cosy'sin (9' — 9) 
— sinl cosI' sin ọ' cos(9'— Q) + sinl cosI sine]. 


En transformant, comme plus haut, cosI et cosT', il vient 
F COSI e) sme |- sin I sin (9'4- 9'— 9) + sinl'sine’ 
Ape NR MEO LCR Pr 1 
+ 2sin Í sin y’ cos (9' — Q) sin? » I' — »sinl'sine'sin*- I |. 
Substituant les valeurs (123), et posant 


1 Es ; E : 
| P't——3i sin (y/4- ?/— 9) sin] + p sin(y/-- #'— 2') sinl 


(128) | 


ovre ; Nat apto 
+ sin (y'+ &/— 9") cos (Q— 9)sinI sin" I— sin (y 4- &/— 9^) sinl'sin*? H 


nous aurons 


(129) PCOS (CE S s.) c r - Y. 


Au moyen de ces valeurs, les équations (119) deviennent 


t > I I I m' 
m (35 a)r [esee #— y — m i| n 


[77] 


P I 

(130) Y,—/fm (5 TÝ 
j kl Æ I T 

Z, — fm (z-a) rir 


et l'on a, pour calculer s', suivant (116), 


* I "x 
(131) $5?— r*$— pr! [cos (+ $'— y — D) + — «| HER 
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Nous allons actuellement procéder au développement de X,, Y,, Z,, sui- 
vant les puissances et produits de «, £, y. 
A cet effet, soit s, la valeur de s’ correspondante à « = o, nous aurons 


(132) sè = r?— 2rr'cos(y/-- &'— y — 9) + 7? 


d’où, en soustrayant de (131), 


puis 


LA | r bd 1 . 3 
Elevant les deux membres de cette équation à la puissance — 5° on 


aura 
(33) 2 1 Jarr' 3.5 BEY sun gru 19.0 CE hjd pk 
Ií ES San l -— eI —— — 7 + 
Ss A LA: 9.4. s d2g,b- xm 065. c : 
d'oü 
MEC. ee Jap o dog MT b Spp 
i BE. ark XOT MN qp ca 7 459 3 
(134) 
d 5.04.9 RS EN 
276618 84 


En mettant ces valeurs dans (130), il viendra 


[X = fm (3 — - gus Mae A. Am (5 — l 


3 arr! 3.5 a2 pp? 3.5. 7 a? r? p^ 
X = j FR S * : br eI eA ... 
+ fm cos (y. + ) á (G së + sj 2.4.6 TUR ) 
(135) 
Spi SiD Wirtpur 
Ta ; /m r je PEG P LO 


£ Sarr. | 3.5ec rr? — 3.5. 3.5.7 e r^r 
—miícat xa 4 WE | 


À I I H I I I 
Y, = fa (25 — pr) sini ex — 9) + m (25 — r8 
0 0 


arp ES atis 3.5 tps 
m' sin (y'4- Ð'— Sls s PACA TT 6L SE ES PT RS PRET © 
(136) + Jin! sin (gl e 70 (Ter Tei = 
1 323 .4w 5. don yu 
ur y a Es : 7 aig —— T 
Nus (5 Lv 2.4 P gr ) 


0 
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f I I I 
Z,—fm'|-5;——Ir-Y 
| 1 J ^ p? 9 | 


(137) M > P 12 - Is 3 5% 
s » x! & wr E. We a k cA m. Ex jh 
E — M s $ > q^ ^r 
+ 2 Cs 2.4 os] 2.5.6 | 33 ) 


Comparons actuellement ces développements avec les formules que 
Wronski a publiées, sous la marque (113), p. LxY. 

En vue de simplifier les résultats, l'auteur prend, pour plan fixe, le plan 
dans lequel se meut la planète m à l'instant considéré, de sorte que l'in- 
clinaison I est nulle, tandis que les quantités I’ et Q’ désignent l'inclinaison 
de l'orbite de m sur celle de m et la longitude du nœud ascendant de m’ sur 
le plan de la méme orbite. D'un autre côté, les longitudes y et y’ des lignes 
fixes étant arbitraires, Wronski fait 


(137 bis) =Y 


Enfin, il néglige dans ses développements les termes du quatrième ordre 
par rapport aux inclinaisons. 
Sous ces conditions, nos formules (124), (126) et (128) se réduisent à 


a—— Sin(y -- ?/'— Q0") sin(y + ^ — 9") sin?l, 
(138) B=— sin(y+ ?'— 9')cos(y + — 9") sin?T', 
(= 2sin(y + $/— 2') sinl. 


Mettant ces valeurs dans les équations précédentes, il vient 


| + 1 1 P: 
| X,— fm (a ji )r'cos/— 9)— 7 
| _\So 1 $9 . 

— = fm! A [7 cos (#'— 9) — r] + ; 4 x 


x r'sin(y + ?/— o')sin(y + — 9')sin?I' —- ..., 


I I : 
3 m (z; — 5) 7” sin (Ð/— 9) 


$5 ] 


I gri. s I I 
C in K sin ('—) sin (y+#—2') + (ss x 7) cos(y-- *—Q) | x 
0 4. z 


x r'sin(y + 9?'— o')sin?I' 4- ..., 


1 1 : å 
Z,— fm! (5 — 5) r'sin(y + &/-—— 2') sinl 
0 


3rr? 


I 1 
M e 


sin*(y + &/— 0') sin(y + P — 9^) sinT'-- .... 


9 
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Ces résultats sont. exactement. conformes avec ceux de Wronski. 

I| est presque superflu de rappeler que les valeurs finales dés X, Y, A 
s'obtiendront en calculant, au moyen des formules (135) à (139), les va- 
leurs fournies par chacune des planétes perturbatrices et faisant les sommes 
respectives. ; 

On a ainsi tous les éléments nécessaires à l'application des méthodes 
d'intégration proposées par l'auteur de la Réforme absolue du savoir humain. 


Il n'est pas sans intérét de faire connaitre que les épreuves de ce travail 
ont été revues, avec soin, par M. O. Callandreau. Nous prions ce jeune et 


zélé astronome d'en agréer tous nos remerciements. xv 


ERRATA. 


Page r7, ligne r, supprimez le mot : mouvement. 
22, ligne 5 en remontant, au lieu de (cos — e), lisez (cosu — e). 
32, ligne 11, au lieu de en négligeant, lisez en considérant comme constante. 
32, lignes 5 et 6 en remontant, mettez dt en diviseur de db. 
(6, ligne 10 en remontant, ajoutez (99 bis) comme marque de l'équation g= y+ «4 — 0. 


PSY ; runt 
(9, ligne 7 en remontant, au lieu de -, lisez - * 
| r a 
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